Prepared for submission to JHEP 



INT-PUB-1 1-057 
CALT 68-2864 



A Formalism for the Systematic Treatment of Rapidity 
Logarithms in Quantum Field Theory 



(N Jui-yu Chiu,"' Ambar Jain,'" Duff Neill,'' Ira Z. Rothstein'" ^ 

'^Department of Physics, Carnegie Mellon University, Pittsburgh, PA 15213, U.S.A. 
^California Institute of Technology, Pasadena CA, 91125, U.S.A. 
E-mail: j ychiuOandrew . emu . edu, ambarSandrew . emu . edu, 
tin dneillSandrew . emu . edu, izrSandrew . emu . edu 

cn 

Abstract: 

Many observables in QCD rely upon the resummation of perturbation theory to retain 
predictive power. Resummation follows after one factorizes the cross section into the rele- 
vant modes. The class of observables which are sensitive to soft recoil effects are particularly 

4^ challenging to factorize and resum since they involve rapidity logarithms. Such observables 

include: transverse momentum distributions at pT much less then the high energy scatter- 
ing scale, jet broadening, exclusive hadroproduction and decay, as well as the Sudakov form 
factor. In this paper we will present a formalism which allows one to factorize and resum 
the perturbative series for such observables in a systematic fashion through the notion of 
a "rapidity renormalization group". That is, a Collin-Soper like equation is realized as a 

psj renormalization group equation, but has a more universal applicability to observables beyond 

the traditional transverse momentum dependent parton distribution functions (TMDPDFs) 

^— I and the Sudakov form factor. This formalism has the feature that it allows one to track the 

^ (non-standard) scheme dependence which is inherent in any scenario where one performs a 

resummation of rapidity divergences. We present a pedagogical introduction to the formalism 
^ by applying it to the well-known massive Sudakov form factor. The formalism is then used to 

study observables of current interest. A factorization theorem for the transverse momentum 
distribution of Higgs production is presented along with the result for the resummed cross 
section at NLL. Our formalism allows one to define gauge invariant TMDPDFs which are 
independent of both the hard scattering amplitude and the soft function, i.e. they are uni- 
versal. We present details of the factorization and resummation of the jet broadening cross 
section including a renormalization in p±^ space. We furthermore show how to regulate and 
renormalize exclusive processes which are plagued by endpoint singularities in such a way as 
to allow for a consistent resummation. 
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1 Introduction and Motivation 

Many observables in high energy collisions suffer from poorly behaved perturbative expansions 
due to the existence of large logarithms. Logarithms of fixed scales, such as masses, are easily 
handled by standard renormalization group procedures within the confines of effective field 
theories. However, when one is interested in less inclusive observables, it is often the case 
that large logarithms of kinematic factors can arise. In particular, when one is interested in 
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studying corners of phase space large hierarchies can induce large logarithms. In such cases 
resummations allow us to maintain control over theoretical errors [1]. 

Within the vast class of observables which require resummation there is a large sub-class 
which are technically more challenging to handle than others. These observables correspond 
generically to those for which the recoil of a collinear jet due to soft radiation is an order 
one effect. We will use the acronym SRSO for Soft Recoil Sensitive Observables. A classic 
example of an SRSO is jet broadening [2, 3]. Another set of SRSO's are transverse momentum 
distributions when pr/Q l, where^ logarithms of the ratio px/Q can invalidate fixed order 
results. 

In this paper we will be employing effective field theory (EFT) techniques to factorize 
and resum large logarithms in these SRSO's. Traditional EFT's allow for resummation in a 
systematic fashion by morphing these large logarithms into logarithms associated with UV 
divergences and then applying standard renormalization group techniques. However, such 
traditional methods where one sums logs of invariant mass scales, are insufficient for SRSO's. 
Indeed, as we will see below for such observables, not all the large logarithms are associated 
with UV divergences in the effective theory. 

Standard SCET methods breakdown when the modal decomposition in the EFT involves 
multiple fields with the same invariant mass scalings (such theories fall under the rubric of 
what is known as SCETn [4]). These cases exactly correspond to SRSO's since the soft 
radiation, in light-cone coordinates, has momentum scaling (A, A, A), where A is the small 
power counting parameter, while the collinears scale as (1,A^,A). For such observables one 
can run into a new type of divergence which is not associated with singular behavior in 
the UV or IR, but with limits of large rapidities, as was elucidated in [5]. These "rapidity 
divergences" have been studied extensively outside the realm of EFT, albeit perhaps with 
differing nomenclature, especially within the context of transverse momentum distribution 
functions and Sudakov form factor [6-8]. 

The purpose of this paper is to present a renormalization group program for the treat- 
ment of rapidity logarithms. In particular, we show how one can regulate and renormalize 
rapidity divergences, and then using a renormalization group technique resum the associated 
logarithms, all in a systematic fashion within SCET. To facilitate this procedure we introduce 
a regulator that necessarily breaks boost invariance in order to distinguish between modes 
which share a mass shell hyperbola. The regulator preserves eikonal exponentiation and man- 
ifest gauge invariance in each sector. However, our formalism is also applicable using other 
regulators which break boost invariance, as will be discussed below. Once one sums over the 
soft and collinear sectors, the rapidity divergences cancel at each order in perturbation the- 
ory. The regulator introduces a new scale, v, which leads to a rapidity renormalization group 
(RRG) flow. The solution to the ensuing differential equation has the effect of summing the 
large rapidity logarithms which existed in the full theory calculation. 

^We will always assume that pt ^ ^qcd, so we are away from the "forward region". 
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Outline of this paper 

We begin in (2) by defining the notion of rapidity logarithms (divergences) and present the 
necessary and sufficient criteria for their existence. In (3) we present a physical argument 
which allows us to isolate the type of generic observable where one expects them to arise. 
In (4) we demonstrate how to regulate the rapidity divergences in context of the massive 
Sudakov form factor with massless external lines. This section also illustrates how one can 
resum the rapidity logarithms using the rapidity renormalization group (RRG). We then 
show how the regulator is applied to generalized soft and collinear jet functions. In (5) we 
apply the formalism to the transverse momentum distribution in Higgs production at small 
Pt, resumming logarithms oi p±/mh to next to leading log order. This section includes the 
definition of a gauge invariant transverse momentum dependent parton distribution function 
(TMDPDF). At the end of this section we compare our results with previous works on the 
subject. Section (6) gives a factorization theorem for jet broadening including a next to 
leading log result for the resummed cross section. Since we have attempted to make this 
section available to readers not interested in the Higgs production section there is some formal 
overlap with the previous section. Finally in (7) we show how our formalism can be utilized 
to renormalize end-point divergences in exclusive decays. We conclude with a summary. 

2 Rapidity Divergences 

We define a rapidity divergence as arising from momentum region where the invariant mass 
k'^ is held fixed but the ratio kj^/k^ (or k^/k^ ) diverges, where k± are light-cone momenta. 
Rapidity divergences are not IR in origin, as they do not show up in the full theory, nor should 
they be thought of as UV divergences since they can arise from either the upper or lower limit 
of an integral, as will be shown later. The existence of the divergence stems from the fact that 
to preserve manifest power counting the EFT must be multipole expanded [9] which leads to 
eikonal propagators. Furthermore, while rapidity divergences arise in factorized IR sectors of 
the theory, i.e. collinear or soft, the sum of EFT sectors will have no rapidity divergences. 
Which is to say that rapidity divergences arise as an artifact of factorization. However, in 
order to resum logarithms in attempting to save perturbation theory, factorization, and thus 
rapidity divergences are inevitable for SRSO's. In this sense, they are entirely analogous to the 
traditional UV divergences which arise in factorization and are necessary for resummations. 

To demonstrate how rapidity divergences arise, let us consider an integral of the following 
form 



which may arise when transverse momentum is measured in the real radiation. Q is the scale 
of hard scattering and /u^ is the relevant low energy scale. Suppose that the n-collinear and 
soft modes have kj^/Q or order 1 and A ~ [ii^jQ^ respectively. Then this integral ranges 
over both mode regions. To factorize the integral into rapidity regions we introduce a set of 
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cut-offs 
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corresponding to soft and collinear contributions respectively. In the effective theory, the 
cut-off should not be finite, to preserve power counting. Or put differently, each sector should 
depend only upon one relevant scale, which follows after performing the multipole expansion. 
Taking the limit, Q » A >> iil, the EFT result reads 



We see that we generate a set of divergences which only cancel in the sum of the sectors. 
It should also be emphasized here that these divergences are not regulated by dimensional 
regularization, a necessary but not sufficient criteria for rapidity divergences. 

Obviously not all observables in QCD will generate rapidity divergences in SCET. A 
necessary condition for their appearance is that the observable under consideration receive 
contributions from modes with parametrically distinct rigidities, but whose invariant mass 
are of the same order. The prototypical observables of this type are transverse momentum 
distributions. When px ^ Q, Pt plays the role of fj^i in the example above and rapidity 
logarithms arise which must be resummed. As such, we will consider in this paper two 
observables, namely Higgs production at small transverse momentum, and Jet broadening. 
Rapidity divergences also occur in exclusive processes. In [5] the authors show that the end 
point singularity problem arises as a consequence of rapidity divergences. We apply our 
formalism to this issue and show how one can systematically renormalize the divergences in 
these exclusive decays. 

3 What Theories Give Rise to Rapidity Divergences: SCETj vs. SCETn 

SCET [10-12] (Soft Collinear Effective Theory) is a formalism designed to separate scales in 
high energy scattering processes for which the hard scattering scale (Q^) is much greater then 
the scale of hadronic physics. Here we will not review SCET but only illuminate the points 
that are germane to the main thrust of this paper. In particular we are interested in the 
critical issue of (not) double counting regions of phase space. That is, how does one cleanly 
separate (factorize) the modes which give rise to IR singularities. 

As in any well defined EFT, the scale separation is made manifest at the level of the 
action which systematizes the power counting. Power corrections can be included by adding 
operators which have definite scalings in powers of hl/Q, where fJ-L some low scale of 
interest. SCET, like its cousin NRQCD, is a "modal" theory whereby fields are decomposed 
into a set of sub-fields each of which has momenta with definite scalings. For instance, in 
SCETi the gluon field is written as 




(2.3) 
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where are collinear fields whose momenta scale as Q(l, A^, A) and Q{X^, 1, A) respec- 

tively. While A^^ is an ultra-soft (US) field whose momentum scales as Q(A^, A'^, A^), where 
the power counting parameter is A = [^l/Q and /^l is the relevant low energy scale^. To 
cleanly separate the scales it is helpful to use a dynamical label formalism [14] also used in 
SCET [11]. We re- write the full QCD field as 

i,ix)= (3.2) 

n-p,p± 

The purpose of this rescaling is to insure that all derivatives acting on ^ scale as A^ . Note that 
(3.2) is written as a sum not an integral. One tessellates the space of possible large momenta 
into bins whose dimensions scale with the size of the residual momentum of order A^ . 

The Lagrangian interactions can change the large momentum components of fields. In 
particular, collinear gluons can split, changing their large light cone momentum. This im- 
plies that there are loops in which one must sum over labels. It is then natural to ask what 
happens in the label sum ^ when one of the labels becomes parametrically small^ and modes 
begin to overlap. The existence of these overlap regions, when not treated properly, obscures 
the physics underlying the effective theory calculations. This is perhaps simplest to see in 
NRQCD, where the existence of the overlaps leads to pinch singularities, as well as the in- 
ability to clearly distinguish between IR and UV divergences. These points were made clear 
in a seminal paper by Manohar and Stewart [5]. In this paper the authors show that one 
may exclude these overlap regions by taking a diagram involving a particular mode, Taylor 
expanding it around the region of the complementary mode one is trying to exclude, and sub- 
tracting this contribution from the original diagram. Doing so eliminates the aforementioned 
problems. It has been shown that in SCETj this "zero bin" subtraction [16-18] is often 
equivalent to dividing by matrix element of Wilson lines in non-SCET perturbative QCD 
factorization formulae [19], which can often be identified with the inverse of the soft-function 
in the factorization theorem. In this paper we will also have occasion to comment on the role 
of the what we term the "soft-bin" subtractions, and draw a distinction between the soft-bin 
and zero-bin, where by zero-bin we always mean an ultra-soft scaling (see appendix (B)). 

In SCETi all divergences can be regulated using dimensional regularization and/or off- 
shellness, and all modes have distinct virtualities (collinear modes in differing directions being 
the exception '''). Zero bin subtractions are relatively simple to utilize to insure that there 
is no double counting. Thus we should expect that rapidity divergences should not be an 

■^The ... represent non-linear terms in the fields that are needed to insure that gauge transformations do 
not mix orders in the power counting [13]. 

^The fact that these are sums and not integrals is a consequence of separating momenta into labels and 
residual momenta. This can be thought of as a grid in which labels give the coordinates of a box (or "bin" ) 
whose size is of order of the residual momentum. One can combine the sum over bins with the residual 
momenta integrals when performing loop calculations. See [15] for a discussion. 

*In SCETii labels between modes can also overlap by becoming large. 

^The soft zero bin will eliminate any such overlap. Given the non-existence of a soft mode in SCETi such 
an overlap is necessarily absent. 
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issue in SCETj , and indeed this is in fact the case. Note this is not to say that integrals of 
the form of (2.2) wih not arise in SCETj . In fact they are ubiquitous, however, they should 
not be interpreted as rapidity divergences, as they will cancel within each sector. That is 
to say, there can, and will be such divergences in the collinear (or US) sector, but when one 
sums over graphs, including zero bins, these divergences will cancel as they must within each 
sector. A classic example arises in the one loop correction to the parton distribution function 
(PDF). In this case in calculating a real correction one encounters an integral of the form 



where z is the momentum fraction carried by the struck parton. The integral diverges at 
the upper end point of integration, where the incoming parton and and the struck parton 
carry the same momentum. To regulate this integral one must introduce a new regulator 
since dimensional regularization is insufficient. Nonetheless the divergence arises when the 
emitted gluon goes ultrasoft, and thus by the Kinoshita, Lee and Nauenberg (KLN) theorem 
we expect it to cancel with the corresponding virtual diagram. Thus the divergence cancels 
within the collinear sector itself. Note that in this case the zero bin did not play a role, since 
they actually cancel in this calculation, again by the KLN theorem. An illustrative example 
of the cancellation of divergences, which one might have thought were rapidity divergences in 
SCETi can be found in [20]. In summary, SCETj rapidity divergences do not pose a problem 
because there is no issue in distinguishing modes of identical virtuality. 

For certain observables, such as transverse momentum distributions with p± <C Q, SCETj 
is not the proper effective theory. The reason is that we must account for real radiation with 
momenta that scale as {p± ,p±,p±)- Given that p± is the IR scale of the theory the collinear 

2 2 

modes scale as either {Q, '^,p±) or {^,Q,p±), and both the soft and collinear mode have 
the same invariant mass. As such, these modes can be interchanged by a boosts and the 
only real distinction between them is in their relative rapidities. In such cases, the EFT is 
called SCETii . The need for SCETn was first pointed out in the context of exclusive B 
decays [4] where it was necessary to introduce a second effective theory below SCETj. The 
equality of invariant masses of the modes in SCETn leads to complications in the factorization 
of physical observables since one must break the boost invariance of operators to cleanly 
distinguish between sectors. As was first pointed out in [21], the process of factorization in 
SCETn can lead to additional divergences in sectors that can not be regulated by dimensional 
regularization, or off-shellness. These divergences will not cancel within each sector as they 
do in SCETj. They will cancel only when we sum over sectors, but the lack of cancellation 
within a sector changes the RG structure of theory. In fact it is the lack of cancellation that 
allows for the resummation of large rapidity logarithms . 

Boost Invariance in SCETi and SCETn 

In SCET we make a convenient choice of frames where relevant modes are (ultra)soft or 
collinear, breaking the full Lorentz symmetry of QCD. However, SCET is invariant under 




(3.3) 
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boosts along the light-cone direction. This residual boost symmetry of SCET is called RPIm 
[22]. For problems relevant to SCETj this symmetry is preserved in each matrix element 
belonging to the collinear or ultrasoft sectors. The natural distinction between the collinear 
and ultrasoft sectors of SCETi comes from their parametrically differing invariant masses^. 
Since dimensional regularization breaks dilatation symmetry, it is sufficient to distinguish 
between collinear and ultrasoft modes^. Since dimensional regularization preserves boost 
invariance and is the only regularization required to separate modes in SCETi, each sector 
individually remains RPIm invariant. 

In SCETii there is no distinction between the invariant masses of the soft and collinear 
modes and they can be interchanged via boosts. Dimensional regularization does not distin- 
guish between these modes, so we should introduce a regulator which accomplishes this goal 
by breaking the boost invariance along the light-cone direction. Due to the soft recoil, "jet 
functions" describing collinear radiation are not exactly aligned with the preferred light-cone 
direction used in factorization, in that they carry transverse momentum w.r.t. the light-cone 
direction^. A jet function will usually depend upon the transverse momentum and the large 
light cone component of the momentum Q± carried by all the collinear particles constituting 
the jet. While the transverse momentum is boost invariant, Q± is not and hence one shouldn't 
expect the jet function describing collinear radiation in a physical process to be boost invari- 
ant. In contrast, the jets in SCETj are aligned with the preferred light-cone direction since 
ultrasofts cannot recoil jets in transverse momentum. 

In a back-to-back jet scenario for SRSO's, like the jet broadening event shape, a small 
boost will reduce the number of particles in one jet and increase the number of particles in the 
other jet while keeping the number density unchanged in the soft region, when averaged over 
all events. Thus it is expected that in a renormalized factorized cross-section boost invariance 
will be broken in each sector and will only be restored when all sectors are added. In problems 
pertaining to SCETj, boosts alone cannot interchange ultrasoft and collinear excitations and 
hence RPIm must be preserved in each sector. We see that in order to factorize one must 
distinguish between soft and collinear radiation, and RPIm must be broken via regularization 
of SCETii matrix elements. This point was emphasized in [24]. 

4 Rapidity Divergences in SCETn 

To understand the nature of rapidity divergences we consider how they arise in the effective 
theory. Thus we will begin by considering perhaps the simplest case where such divergences 
arise, the Sudakov form factor. In particular we will renormalize and resum the logarithms 

^Distinction between two collinear sectors arise essentially from different light-cone directions, so they need 
not have different invariant mass for distinction. 

'^An ultrasoft mode can be transformed into a collinear mode by a dilatation and a boost. 

*For example, jet functions in the jet broadening factorization theorem [23] carry non-zero transverse 
momentum w.r.t. the thrust axis. 
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in the on-shell, massive gauge boson, space-like Sudakov form factor [19] ^ ^o. As opposed 
to the massless form factor, which in isolation, is unphysical, the massive case is IR safe and 
in principle observable. This form factor is also relevant to for summing large electro-weak 
corrections at energy far above the gauge boson masses. 

We defined our power counting parameter as A = M/Q. The factorization formula 
is composed of decomposition in terms of modes which can contribute to the non-analytic 
structure of the matrix element. The relevant modes for this observable are the soft with 
light-cone momenta scaling as (A, A, A), and collinear and anti-collinear with momenta scal- 
ing as (1,A^,A) and (A^,1,A) respectively. There are no ultra-soft contributions since their 
momentum (A'^, A^, A^) will decouple from all the other lines. Given these modes, this defines 
an SCETii process. The jets will recoil against the soft virtual emissions and hence, despite 
the exclusive nature of this process, it is still an SRSO. However, in the Breit frame the net 
transverse momentum exchanged between the jets must vanish. 

The factorization of the massless Sudakov form factor in SCET was performed in [26]^^. 
At leading power we have 

= u{pn)lj^u{pn) F{Q\ M2) « {pn\LWnShj^C{n ■r,n- V)SnWl^n\Pn) (4.1) 

= H{Q^, fi)Jn{M- /X, u/Q)-f^Jn{M- /x, iy/Q)S{M; /i, u/M) . 

In the last line, we have factorized the form factor F in terms of SCETn matrix elements 
Jn, Jn,S which are defined as 

S{M;jj,ij/M) = (0 I S^S I 0) 

J„(M;/i,Z./Q) = {pn\^nWn\0) 

JniM;f,,u/Q) = {0\WlCn I Pn) (4.2) 

S and W are Wilson lines composed of soft and collinear lines respectively. 

Finally we must mention the cumbersome issue of the Glauber modes. These modes, 
which scale as (A^,A^,A) will in general contribute IR singular pieces at the level of the 
amplitude. While they have been shown to cancel in certain processes outside the realm of 
effective field theory [27, 28] , a systematic treatment of such modes within a self-consistent 
EFT treatment is still lacking. Here we will assume, as do all SCET treatments, that Glaubers 
will not contribute. 

®The logarithms in the Sudakov form factor are distinguished from "Sudakov logarithms" which can arise 
in running currents in that they contain rapidity divergences. 

The Sudakov form factor can be regulated by going off-shell in which case the systematics may change 

[25]. 

^^This factorization was for the massless case where gauge invariance uniquely fixes the form of the factor- 
ization theorem. However, the gauge boson mass does not alter the result. Note that this result was formal in 
that the IR scales were not clearly delineated. As such, it was not sharply defined to be living in SCETi or 
SCETiiand the factorization formula contained both soft and ultra-soft Wilson lines. 

^^For a discussion of these modes within the context of SCET see [26, 29, 30]. 
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Figure 1. The mass-shell hyperbolae showing the distinction between the different sectors [5]. The 
separation between soft and coUinear modes is arbitrary and leads to rapidity divergences. The soft 
sector has two distinct rapidity (UV) divergences that must cancel with rapidity (IR) divergences 
arising from the collinear sector. 




Let us now see how factorization of the soft from collinear modes leads to rapidity diver- 
gences. Consider the full theory one loop vertex correction. The relevant scalar integral is 
given by 

/" n 1 1 1 

If J [d k] _ j^j2^ l^j^2 _ ^ . f^f^ . _^ (^]^2 _ fi . j^j^ . _ ^g-) 

This integral is finite in UV as well as the IR. In the effective theory there are three 
contributions. A soft integral coming from taking the limit — )• {M, M, M) 

Is = j K/c] _ ^2) A: + ie) {-n ■ k + ie) ^^'^^ 
and two collinear integrals (/n, In) of the form 

- M2) {k'^ -n-kn-pi+ie) {-n ■ k + ie)' ^^'^^ 

Given that the full theory graph is IR finite, so must be the sum of the effective theory 
graphs. Let us consider the soft graph integrating over k±. 

Is ^ f[cPk]{n-kn-k-M'^)-'^' ^ ^ 



[—n - k + ie) {—n - k + ie) 

(4.6) 

We see that the relevant region of phase space lives on the hyperbola n ■ kn • k ^ , shown 
in figure 1. Off the hyperbola the integral becomes scaleless. Given this restriction, we note 
that the integral diverges when the rapidity (n ■ k/n ■ k) approaches infinity or zero. These 
divergences are not regulated by dimensional regularization and correspond to the rapidity 
divergences that arise when the soft integral overlaps with the two collinear rapidity regions. 
This is illustrated in figure (4). On the other hand, if we consider the collinear n diagram 



- 9 - 



we see that it only has divergences associated with the hmit where (n ■ k/n ■ k) approaches 
infinity, and similarly with (n — )• n) for the In collinear integral, since there is only one border 
between a collinear sector and the neighboring soft sector. 

There are multiple ways of regulating these rapidity divergences. One can go off the 
light cone by setting / [7, 19], use an analytic regulator [31], or a "delta" regulator as 
was done in [32]. Choosing a regulator determines how much of an overlap there is between 
modes. For instance, in the case of a delta regulator, where one shifts the eikonal propagators 



one must perform a soft-bin subtraction to generate the correct result in the effective theory. 
In fact, the authors of [32] showed that the sum of the integrands, once properly soft-bin 
subtracted leaves a finite integral with no rapidity divergences. With an analytic regulator 
the soft function vanishes explicitly. In this case there is no double counting as half of the 
soft contribution comes from each of the collinear sectors, and thus there is no zero bin. 

While physically it seems clear that a sensible rapidity regulator should cancel in the 
sum over sectors, we should have a proof of this assertion. A direct proof follows noting that 
if all of the regulated EFT diagrams arise from an asymptotic expansion of the full theory 
diagrams, then given that the full theory has no rapidity divergence the finiteness of the EFT 
sum then follows. By this reasoning the delta regulator must also cancel in the sum over 
sectors, as will any rapidity regulator if we assume that it preserves the equality between 
the full theory integrals and their asymptotic expansion in regions. It is important however, 
to recall that the method of regions is distinct from EFT in that, in the latter, it is not 
necessarily true that there is a one to one correspondence between the IR of a full theory 
diagram and a corresponding diagram in the EFT. All that is necessary is that the sum of 
the EFT graphs reproduces the IR of the full theory. Such cases arise when one uses the 
equations of motion in the effective theory to remove redundant operators. 

4.1 Regulating the Rapidity Divergences 

Given that EFT's are created to sum logarithms we would like to be able to regulate the 
theory in a way that makes an RG treatment manifest. There are multiple ways in which 
to regulate the rapidity divergences, and the formalism developed here can be applied using 
any sensible choice, such as the delta regulator [32]. Here we will concentrate on the regu- 
larization introduced in [23], where we utilized a rapidity regulator which is closely related 
to dimensional regularization. It is implemented by modifying the momentum space Wilson 
lines in the following fashion. 



1 



1 



(4.7) 



n • k 



n - k + A 




,2 




perms 




2r, 



rj/2 



(4.9) 



perms 
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We have introduced a new dimensionful parameter v which will play the role of an effective 
rapidity cut-off. Here is the momentum operator and we have essentially regulated the 
longitudinal momenta, and since I2P3I — )• \n ■V\ in the collinear limit. Note the differing 
powers of r/ in the soft and collinear Wilson lines. The appropriate power is fixed by ensuring 
that the rapidity divergences cancel to all orders which we shall show below. Alternatively, 
and equivalently, the power is fixed by regulating the full theory diagram and taking limits of 
the integrand. The relative factor of two comes from that fact that for a given gluon line in 
the full theory there are two soft eikonal vertices (connecting the two eikonal lines) relative 
to the one collinear eikonal vertex. We have also introduced the bookkeeping parameter w 
for convenience, which eventually will be set to one. It will play a role when we derive RG 
equations. The g subscript on the momentum (label) operator will only play a role when we 
consider going to higher orders as is explained in section (4.4) and appendix (A). 

With this regulator the effective theory will have divergences in both the r\ and e go 
to zero limits. The order of the limits is crucial to sensibly renormalize the theory. Given 
our physical arguments regarding the nature of the rapidity divergences, the proper order of 
limits must be: r/ — )• 0, then e — )• with r//e" — )• for all n > 0. The physical reason for 
this ordering is clear since we must remain on the invariant mass hyperbola when we take the 
rapidity cut-off to its limit. To see how this works in practice let us evaluate the integrals Is 
and In using this regulator. 

The Is integral is most simply evaluated by first doing the /cq integral by contours. The 
result, after repristinating the expression with the coupling, group theory factor and the 
relevant numerator for the Sudakov form factor, in Feynman gauge, is given by 



Is = -g^CF{e^^'2 



-r,-2 -5/2X (J^Y' (J^y r(l/2 - 7?/2)r(e + 7?/2) 



Expanding first in r/ and then in e we find 

,2e 



Is = g^Cp 



Mj 



+ 



1 fH^) , .2. ^ 



47r2 



Similarly, the collinear integral In is given by 



In = g^Cp 



n ' \ ^ ln( — ) ln( ) + ln(-^) H 1 + ln( 

Sn^rj 47r2 V ^ n ■ pi' ^M' 2e\ n ■ pi 



(4.10) 



1 

96 
(411) 



+ 



(4.12) 



1 

48 



and In by replacing n ■ pi with n ■ p2- Summing the sectors we find 



Is+In+In = g Cf 



^\^ + — + - +ln2(-^ +2 n( — ) + 21n — In — + 1 

47r2 \^2e2 e e ^M' ^M' ^ M j 96 

(4.13") 



w has been set to one, and is utilized below when we derive the renormalization group equation. We have 
also absorbed the MS factor into fj, to simplify the expressions. 
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where we have used n-pi = n-p2 = Q. We see that the r] (rapidity) divergences vanish, there 
is no dependence on the scale and the answer is boost [22] invariant. 

In addition, note that the soft-bins are all scaleless and vanish. However, this does not 
mean that they should be ignored, as explained in appendix (B). Indeed, as emphasized in [5], 
these subtractions can play a crucial role in being able to discern IR and UV singularities. In 
the case of the 77 regulator this scaleless, vanishing, soft-bin contribution has the effect shifting 
the rapidity cut-off to its proper place. That is, if we are regulating the a collinear integral 
the effect of the soft-bin will (formally) to shift the cut-off to its proper place separating the 
collinear from the soft. 

4.2 The Rapidity Renormalization Group 

An advantage of the regulator we have introduced is that it allows one to write down a 
renormalization group equation in a rather straightforward manner. We begin by examining 
the Sudakov form factor of the space-like current in terms of the SCETn fields, 

Jf, = HiQ^fi)Jn{M; ^, iy/Q)-f^JniM; /i, u/Q)S{M; /x, i^/M) (4.14) 

The one loop values of matrix elements Jn,Jn,S defined in (4.2), are given by (4.11,4.12). 
The renormalizaton group follows from the set of equations 

(J„, 5)^^-^^ = -j^iJn, S)""^'" = 0. (4.15) 



filn[;u] (iln[i/ 

Moreover the independence of and leads to 

J d d 
dln[fj] ' d ln[i/] 



0, (4.16) 



which is of course true for any observable not just the Sudakov form factor. 

Defining the anomalous dimension under and v variations as (7/^,7^) respectively, such 

that 

-C'--Z-l,^^Z„^s, (4.18) 
equation (4.16) imposes the constraint 

which holds for any observable of interest. Z is an integer whose value depends upon whether 
we are considering an amplitude or the square of an amplitude. For the Sudakov form factor 
Z is either 1 or 2 (see below). The last equality comes from the consistency of ^u-anomalous 



- 12 - 



dimension with the hard anomalous dimension which is hnear in the logarithm with coefficient 
Tcusp- The universal relation between the collinear c and soft S anomalous dimension 

- 2Zc = Zs (4.20) 

follows automatically from the v independence of the hard function, as will be discussed 
below. 

Let us now apply the RRG to the Sudakov case we studied above. Since our regulator 
allows us to define the jet and soft functions independently we may renormalize them in 
standard fashion by absorbing ^ and ^ divergences in the renormalization constants, and 
then run renormalized quantities individually. We define the renormalization factor Z„, Zs 
via 



J, 



R 



yl/2y-l jB qR 

^ij) -Jn 



-1 cB 



(4.21) 

where corresponds to bare quantities and to renormalized. Then using our result from 
above, at one loop we have 

v2e 



Zs = l 



Zn 



1 + 



47r2 
47r2 



2e2 



In^ 

V 



^''^m (^) 

2r] 



2e 



1 ^ . V 

+ — 1 + In 

2e \ n • pi 



where Z^ is wave function renormalization which is the same as in full QCD. 



The [1 anomalous dimensions are given by 



g(/^)^gF 

167r2e 



(4.22) 



(4.23) 



4^2 



In 



7^ 



47r2 



3 , V 

- +ln^ 

4 n- pi 

3 ^ u 

- +ln 

4 n- p2 



(4.24) 



47r2 



In^ 



/2 • 



As a consistency check see that 



+ + 7^ 



47r2 V 2 



(4.25) 



where 7// is the anomalous dimension of the hard matching coefficient. 

The calculation of the v anomalous dimensions necessitates care. The bare book keeping 
parameter is v independent and thus, in analogy with the coupling the "renormalized"^^ 
w obeys 

(4.26) 



d 

v—w 
ov 



r] 



-w 



It is important to remember that w is not a coupling, but strictly a calculational tool. 
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we find at one loop 




87r2 M2 ' 
47r2 M2 



(4.27) 



These correctly obey the consistency equation 



(4.28) 



Both the large logarithms, due to large invariant mass ratio and large rapidity ratio, can 
be resummed by the RG equations 



The relation (4.19) guarantees that the /i and v evolutions commute, hence, the evolution 
in ^-v plane is path independent. However, care must be taken when solving the i^-RG 
equation, "fy contains terms of form a" (/i) In™ (|u/M) with m < n. For instance, one can 
see from Fig. 2 that the one loop result will be multiplied by a series of logarithms of the 
form ^,„[/3oas ln(/x/M)]". These logarithms can be large if ^ » Af, for example, and would 
require resummation. This is easily obtained by solving the consistency relation (4.19) up to 
the required order in perturbation theory. 



where integration constant is fixed by the fixed order calculation of anomalous dimension 
and corresponds to its non-cusp piece. From eqn. (4.27) we see that non-cusp piece is zero 
at one loop. Eqn. (4.30) completely fixes the logarithmic (/i) structure of "fy to all orders 
in perturbation theory when expanded in as{^). If we had calculated "fy to higher orders 
we would see these logarithms explicitly. Thus, it constitutes a check on the higher order 
calculations. In its integrated form resums the set of diagrams which renormalize the 
coupling, which in the Abelian case, arise from the bubble chain shown in figure 2, thus 
taking into account the running of a^. Fixed order form of suffices when evolution is done 
along path 1 shown in figure 3 with /ij ~ ~ M <^ fif Vf. However, the integrated 
form (4.30) is required when evolution is done along path 2. Since fif ^ M there are large 
logarithms in 'j^, that require resummation in addition to the rapidity logs. In figure 3, U 
and V are the evolution factors in /j, and u respectively and /ij, fj are the scales for the initial 
conditions. The notation U {fj. f , fii; Ua) implies running fj, from fii to fij at fixed ly = Ua] 



U^{Jn,S)=j:''{Jn,S). 



(4.29) 




(4.30) 
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Figure 2. Coupling renormalization (Abelian) contributes to 71, and is missed in the fixed order one 
loop result. 



y(z//, z/^;/x/)^fi nal 



[pat 



path 2 



pa th 1 f 



initial V{vf, vi] Hi) 



H >l> 



Figure 3. Two alternate paths are shown for evolution in fi-i^ plane. Due to independence of and 
v scales evolution along cither path will yield the same result. 



similarly for V{vf, h'f, ^a)- Along path 1, we have chosen to run first in u and then in /i. Path 
2 shows the alternate choice and should yield the same result, thus 



(4.31) 



To ensure this in practice, we must use the resummed form of when calculating V{vf, Ui] jif). 

Notice that these anomalous dimensions depend upon the "low" energy parameter, M, 
which normally would, and should, not show up in the expression for an anomalous dimension. 
However, we must recall as far as the rapidity divergences are concerned M is not a low 
energy parameter, but just the invariant mass of the hyperbola along which the rapidity 
renormalization group flows. 

To sum the large logarithms we first identify the natural scales for the Hard, Soft and 
Jet Function which are given by (/U^), (/^s, J^s) and {fij,i>j) respectively. Numerically they 
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can be read off from (4.11,4.12) 



(4.32) 



To eliminate the large logarithms we may run in both /i and u to some fixed scale, while 
evaluating the fixed order functions at their natural scales. That is, we may write 



S{^i,v) = Vs{i^,i^s;fJ')iUsifJ;fJ's;'^s)S{fj.s,t^s)) 



(4.33) 



where Un,s and Vn,s are respectively /i and v evolution factors for jet and soft functions. In 
(4.33) we have chosen to run first in and then in u. We could equally well have switched 
the order leading to the same result. Note that in the ordering of eqn. (4.33) we are required 
to use the integrated form of of eqn. (4.30) in order to resum all the large logs due to the 
running coupling. We get. 



Us{^J',^^s■,l's) = exp 
ys{v,vs\l^) = exp 

Uj{n,fij;uj) = exp 
Vj{u,uj;fi) = exp 
UH{tJ',fJ-H) = exp 



SttCf 



1 



1 



In 



f^l a{ns) a{us) a(/is) 



/3o 



In 



a{M) J 



2Cf /3 I, fi^j 
Cf , / a(/x) \ 



In 



a(/ij) 



SttCf 



1 



In 



/3o \a{fJ-H) a(/i) a{Q) aifJ-H) 



with 



S{fJ-s,'^s) = 1 + 

Jn{l^J,Vj) = 1 + 



a{lJ's)CF 
vr 

a{fj,j)CF 



ln2(^)-21n(^)ln(^: 



vr 
24 



ln(^)ln(^ 
M n ■ pi 



4 ^M' 



TT^ 1 

12 ^ 2 



(4.34) 
(4.35) 
(4.36) 
(4.37) 
(4.38) 

(4.39) 
(4.40) 



Using relations (4.34) to (4.37) we can explicitly verify the commutation relation (4.31) 
at the order we are working. Equations (4.33) to (4.39) give the resummation for the most 
general choice of scales fi and i^. However, in order to resum all the logarithms, the most 
convenient choice of scales is /U = /ij = ^5 ~ M and v = vj ^ Q. Running with this choice 
of scales only requires running the hard function in /x and soft function in u to the natural 
scales of the jet function. This strategy is shown in figure 4. With this strategy, it is not 
required to use the integrated form (4.30) and the fixed order form of suffices. 

The physics of the RRG flow can be understood from figure 5. A change in the scale v 
corresponds to a flow between the soft and collinear regions. The natural scale for the soft 
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M-- 



soft 



jet 



M 



Figure 4. Simplest running strategy to resum all the large logarithms in the Sudakov Form Factor. 



function is n-/c~n-fc~M whereas the collinear functions sit at the scale Q. To sum the 
logarithms we may slide the cut-off(s) of the soft function up the hyperbola, such that the 
scale v minimizes the logarithms in the collinear sectors. 

4.3 The Necessity for RRG 

The RRG is critical in establishing the /i independence of the resummed form factor. To 
illustrate this we can combine the evolution factors, and present the completely resummed 
form factor as 



F{Q^,M^) = E{fi, fiH, MJ, ^5; i^J, iys)H{Q^fijj)Jn{f^j, uj; M; Q)Jn{fij, Kj; M; Q)S{iJs, Ks; M) 

(4.41) 



where we have made clear that the jet and soft function depend explicitly on the scales M 
and Q. The ji dependence in E is always sub-leading, and would cancel in the exact result. 
We will keep 11 to be arbitrary to show how the /i dependence cancels to any given working 




n-coll. 



Figure 5. Running in v corresponds to flow along the mass-shell hyperbola. 
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order. Minimizing the logarithms in ah sectors can be accompUshed by the choices: 



m = Q f^j = ns = M ^ n (4.42) 

iyj = Q us = M. (4.43) 



Then the total evolution factor then at one loop is: 



E = exp 



8.c,f 1 1 1 ,J-(4]]^AJ^],JQ' 



/32 \a{Q) a{ii) a{Q) \a{Q) J J \a{M) J 



RRG 



(4.44) 



To the logarithmic order we are working, the exponent of the form factor is //-independent. 
We achieved this critically important feature because of the rapidity renormalization group. 
The variation of the hard double logarithms must be canceled by the variation of the jet and 
soft sectors. Given the two scale nature of our SCETn problem, it is not possible to cancel 
the variation of the hard double log without the presence of a large logarithm in the low-scale 
matrix elements. 

Schematically we can write for a generic soft-coUinear factorization of a physical process 
featuring large double logarithms: 

(Tresum = exp [r[a]L^ - 2r[a]L L + ...] /(L) (4.45) 
L = In(^) (4.46) 
L = ln(^) (4.47) 

Where M is the infra-red scale, Q is the hard scale, and is assumed to be of order M (thus 
L:^ L). The function / is the low scale matrix elements of the collinear and soft sector, and 
should have no large logarithms. We have neglected the running in a, but its inclusion does 
not change the substance of our argument. Varying n we find that the leading fi variation in 
the logarithmic power counting (i.e., terms that scale as L6L) is given by: 

ScTresum = crresumS(T[a]L^ - 2r[a]L + ... (4.48) 
= + ... (4.49) 

The exponent has the required form to cancel the leading variation, since 5L = 5L. In general 
terms in the resummed exponent such as r„_ia"L^ get canceled by an RRG resummation of 
the form 2rn-ia'^LL. There are further sub-leading variations that scale as a"L or a"^, but 
these variations are canceled by the matrix elements as the low scale since they involve no 
large logarithm. 
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This form of the exponent is found in both SCETi and SCETn as the hard double loga- 
rithmic terms appear generically when there are collinear and soft modes, irrespective of the 
scaling of the soft modes themselves. SCETj differs from SCETn in that the corresponding 
LL term in SCETj comes solely as a consequence of traditional RG techniques. One can 
see that the low/high-scale mixed double logarithms can be generated through traditional re- 
summations since the virtuality of the various modes are separated by an equivalent amount, 

q2 q2 

namely ~ -rS^ ~ A . Hence in the running from the jet to the ultra-soft sectors, a double 



log of the form LL is generated . Since SCETn is bereft of the invariant mass separation in 
the infra-red physics, one can conclude that something like the RRG must exist in all SCETn 
factorizations to generate the required LL term in the exponent. 

It is important to note that the inclusion of the rapidity logs in the exponent, necessary 
for insuring ^ independence, also leads to an ambiguity in the scale of the rapidity logs. In 
previous SCET treatments, ji independence is achieved by explicitly exponentiating the large 
logarithm found in the low scale matrix elements [33]. However, doing things in this 
way obscures the errors at higher orders in perturbation theory, due to an ambiguity in the 
choice of low scale in the exponentiated single logarithm, /x independence only tells us that 
we must have a ln(/i/A) in the exponent, but the scale A is arbitrary. The dependence on 
A should of course cancel between the matrix element of the exponent (at any given order), 
which is achieved automatically in the RRG. One can exponentiate the rapidity logs by hand, 
without recourse to the RRG, but then it becomes difficult to track the scheme dependence 
(i.e. how varying the choice of A affects the sub-leading pieces). In contrast, RRG provides an 
independent scale v to vary so that size of all the sub-leading logarithms is properly captured. 
We propose to quantify error estimates by independently varying in /i and u va. a suitable 
range and then adding the errors in quadrature. 

Finally, in earlier, non-EFT, treatments of the form factor, the /i independence was 
achieved via the Collins-Soper equation [7]. In the CSS approach to the form factor the 
resummed exponent is written as [19] 



By running in Q, this effectively mixes rapidity and running since the hard function depends 
upon the scale Q. By introducing the scale v we avoid this issue, which allows us to cleanly 
separate the rapidity logs from the invariant mass logs. However, note that this methodology 
will lead to two independent integration constants just as in our formalism. Having two such 
integration constants allows one to systematically control scale variation errors. 

^^This statement is of course contingent on the arbitrary choice of which functions to run. Different choices 
would have these logs appear in other evolution factors. But their appearance is intimately tied to the invariant 
mass scale separation of the effective theory sectors. 

In this work it was shown that there is at most one log in the low scale matching. 



H 




(4.50) 
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4.4 Gauge Invariance and the Structure of Rapidity Divergences at Higher Or- 
ders 

In any approach to renormalization, one does not want the procedure to violate gauge in- 
variance. That the total contribution, i.e. the sum of the soft and coUinear pieces, is gauge 
invariant follows from the same argument used to show the rapidity divergences cancel in the 
effective theory. One first introduces the regulator in the full theory, where the regulator is 
not needed to make integrals well defined. Since the rapidity regulator is inherited by each 
sector (in appendix A, we show how this is specifically accomplished for the ??-regulator) from 
the full theory, it is guaranteed to cancel in the sum of sectors. Gauge invariance follows 
similarly. That is, the full theory is gauge invariant, and under the assumption that the 
effective theory is properly reproducing the infrared physics, then the sum of the effective 
theory diagrams must also be gauge invariant. 

Given the intimate connection between Lorentz invariance the gauge symmetry is might 
seem surprising that the regulated sectors are themselves gauge invariant. Nonetheless, as we 
prove in appendix (A) the regulated sectors are indeed gauge invariant to all orders, in covari- 
ant gauges. This result follows once one introduces the notion of non-abelian exponentiation 
[34, 35], (see [36] for a nice modern approach) which strongly constrains the structure of the 
rapidity divergences. The diagrammatic expansion of any generalized soft function can be 
rewritten as the exponential of a distinct subset of diagrams contributing to the series. That 
is, a generalized soft function can be written as the exponential of two eikonal line irreducible 
graphs with a particular color weight. The sum of these graphs is known as a CWEB. Since 
only a single rapidity divergence can appear in the logarithm of a generalized soft function 
[37, 38], it follows that only a single rapidity divergence can appear in a CWEB, regardless 
of the number of loops involved. Hence CWEBs are minimally divergent with respect to 
rapidity divergences Indeed, this fact is critical in establishing the gauge independence 
of the rapidity anomalous dimension. Given this marginal divergence, the gauge dependent 
piece of the covariant gauge polarization tensor, which is proportional to kfj,k^, will lead to 
an integral with no rapidity divergence. 

To see this, consider an n loop CWEB. Since the marginal nature of the divergence implies 
that there are no sub-divergences, we may perform the n — 1 loop integrals leaving one loop 
integral left over with a gluon attaching to an eikonal line. Given that we have performed 
all the loop integrals except one, the gauge dependent piece of the polarization of this gluon 
must be proportional to the momentum carried by the Wilson line itself. Contraction with the 
eikonal vertex will then cancel the denominator which is the cause of the rapidity divergence. 

In appendix (A) we give details of this argument, and show how one needs to generalize 
the Wilson line regulator at higher orders. This appendix also contains a discussion of the 

^'^For rapidity divergences, this statement is expressing the fact that the anomalous dimensions is at most 
hnear in logarithm associated with the rapidity divergence. There are of course other UV renormalization point 
dependent logarithms in the rapidity anomalous dimension, since a CWEB at higher orders have multiple sub- 
loops. But these logs are predicted by the UV divergences of QCD: see (4.30). 
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delta regulator. In either case, organizing the calculation in terms of CWEBs greatly reduces 
the amount of work necessary. 



5 Transverse Momentum Spectrum in Higgs Production 

The differential cross-section for producing the Higgs boson with fixed momentum in hadronic 
collisions is an observable of obvious relevance. It has been shown that the p± distribution can 
be used as a smoking gun for new physics [39]. Thus having a reliable theoretical prediction 
is a worthwhile enterprise. In the limit where p± <^ rrih this cross section provides another 
important instance where factorization proceeds in SCETn, involving rapidity divergences. 
The Higgs p± spectrum is kinematically and formally very similar to the Drell-Yan p± spec- 
trum. Much effort has been devoted to both understanding factorization in SCET framework 
[40-43] and resummation of logarithms of the form ln(p^/m/i), usually in the context of the 
CSS resummation formalism [44-49]. 

To see why this observable fits into SCETn, let us consider the kinematics. We impose the 
kinematical constraint that the transverse momentum of the Higgs relative to the colliding 
beam be small compared to the Higgs mass, A ~ — <C 1. Thus all final state radiation 
recoiling against the Higgs must fulfill this same condition. Taking the Higgs momenta to 
scale as 



it is simple to see that the on-shell radiation that can recoils against the Higgs scales as: 



where {pc,Pc,Ps) stand for collinear, anti-collinear and soft momentum respectively. We have 
allowed for the possibility of jets in the central region with momentum pcj- These jets each 
have large transverse momenta, but their net transverse momentum must scale as A, and they 
impart little transverse momentum to the Higgs itself. At fixed order in QCD, these jets do 
not appear until NLO in the p± spectrum (or NNLO in total Higgs production). In what 
follows we will prove a factorization theorem that robustly accounts for all such radiation, 
and calculate the resummation to NLL. 

5.1 QCD cross-section 

Incorporating the most recent bounds from the LHC [50, 51] we will assume the Higgs is 
sufficiently light that its dominant production mechanism is gluon fusion. Given this as- 
sumption, we may work within the Higgs effective theory where the top quark is integrated 



Ph ~ mh{l,l,X) 



(5.1) 




Pc ~ mh{l,X^,X) , 
Pc ~ mh{X^, 1, A) , 
Ps ~ m/j(A, A, A) , 
PCJ ~ rnh{l, 1, 1), 



(5.2) 
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out, generating the dimension six operator 



n{x) = h{x)Tv[G^"'{x)G^u{x)]. (5.3) 
The matching coefficient for this operator is known to two loops and is given by [52, 53] 

The differential cross-section in Higgs boson transverse momentum {p±) and rapidity (y) is 
given by 

"^'^ - / d^x ^ {pnPnl nx) 5(y -\^^^) - \A±\'')'HmPnPn) ■ (5.5) 

spins 



dp\dy 8v^S 



Vh is the momentum operator that picks out the Higgs momentum, consider it a derivative. 
\PnPn) is the incoming proton state with momenta Pn and Pn- v is electro- weak symmetry 
breaking scale, and ^/S is center of mass energy. When there exists central jets the cross 
section will become sensitive to higher dimensional operators, but as we will see below, this 
region of phase space is power suppressed. We can simplify (5.5) by writing 

{0\hix)6(y - I In 5(pi - \Vh±\')hmO) 



2 



^'^'^ {27r)5+{pl - ml)e-^^^-6(^y " ^ In ^) S{pl - \p,^\% (5.6) 



the cross-section then becomes 

da _ Cf f d^ph + 2 2x/ 1, Ph' 



dpldy - S.^S I - - I I) ^^-^ - 

J d^xe-'P'^-'^ {VnPnmG''%x)G^,{x)\T,[G''^{Qi)G^p{mVnPn) , (5.7) 



spins 

where 5~^{pf^ — m\) = 0{p^)5{p\ 

5.2 Factorization in SCETn 

5.2.1 Central Jets are Power Suppressed 

To match to the effective theory, one should perform an OPE that matches the full theory 
operator in (5.7) to a product of effective theory operators at the hard scale, as is done in the 
case of inclusive Drell-Yan [54] . As long as one can show that the contributions due to hard 
colored particles crossing the cut in a given full theory diagram is power-suppressed in the ^ 
expansion, then effectively there is no OPE since the currents are still separated by distance 
scales large compared to rrih. If there are no hard partons crossing the cut, then matching at 
the hard scale reduces to matching the full theory operator Tr[GG](x) onto effective theory 
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Figure 6. A contribution to the production cross section with central jets. 



currents 13^"iS'ij_^{x) , where B^"^ wih be defined below. In what follows, we will establish 
that the central jets are power suppressed, and so one can simply match currents. 

At leading order in the p±/mh expansion, the only full theory diagrams that contribute 
must come with a 6{p\) or the ^ singularity, as both of these distributions are leading order 
in the power counting. Delta function and power like singularities are associated with virtual 
and coUinear contributions respectively, neither of which can arise from central jets. When 
matching at the high scale we expand in powers oi p±. Since the central jets only depend on 
the hard scales, one can set transverse momentum to zero in diagrams containing the central 
jets. Thus the part of the full QCD amplitude-squared that produce any modification of the 
Higgs transverse momentum by central jets is power suppressed. 

More formally (pedantically), consider a full theory diagram that contains central jets, 
such as in Figure 6. The diagram has the form: 



Icj 



(2vr) 



;;i^^{kf)j7r^^^{k2)77r^^^iPh-^h)^^'^\Pn+Pn 



(2vr) 



(5. 



The hard contribution to the matching from this diagrams can be obtained by simply con- 
sidering the ki and k2 momenta to be large, and asymptotically expanding the integrand 
accordingly (along with the power counting of the external momenta Pn,Pn,Ph)- Then the 
above integral becomes: 



(5.9) 
(5.10) 



Icj\hard = / d9.6{n.pn - n.ki - n.k2 - n.ph)6{n.pn - n.ki - n.k2 - n.ph) 



dn 



iPH-m,) 



2± ) Pasym {Pn, Pn, h, h, Ph) 



Since all propagators in the integrand Pasym have a hard scaling and the momentum conser- 
vation delta function lacks the scale p±, the integral is completely independent p±. Hence 
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it has no contribution that scales as -y-. Thus the hard contribution from this diagram is 

Pi ^ 

power suppressed. The argument easily generalizes to more complicated diagrams. 
5.2.2 Matching 

Having established that central jets are power-suppressed, we first match the non-local oper- 
ator T-l{x)'H{0) onto the product of SCET currents: 



(M|Tr[G^^(x)G^,(x)]Tr[G°''(0)G„;3(0)]|M) 

= Y,H{mh){pnPn\ (BTA^)B-n±^{x))(Btl{0)B'^^M)\PnPn) + 0{X) . (5.11) 



The interaction is non-local only along the light cone and the transverse directions as there 
are no hard partons crossing the cut^*^. The hard matching can depend only upon ruh and is 
given by 

H{mh)=4mi\C{mh)\\ (5.12) 
where C{mh) is the matching coefficient from current matching^^. Furthermore 



= -Tr 



Wl{x)iD>'^^Wn{2 



and 



n ■ An{nX)dX I . 



Wn is a collinear Wilson line in the fundamental representation defined in x-space by 

Wn{x) = Pexp (^j 

S^' {x) is a soft Wilson line in the adjoint representation 

Sf{x) = Pex.pU n- As{n\)d\ 

\J —oo 

Now we factorize the matrix element: 

(5.11) = F(m,)(p„|i?i (x)i?^'j:(0)|p„)(p^|i?^'j:(x)i?^'l(0)|p^) 
(0|Sr'(x)5f (x)5^^'(0)5f (0)|0) 

= (i^2 -'l)2 (Pn|tr[<^(x)ij-^(0)]|p„)(pdtr[i?^^(x)i3;;^(0)]|p, 

{0\Snx)Sf{x)St'{0)S^\0m 



(5.13) 
(5.14) 

(5.15) 



(5.16) 



(5.17) 



(5.18) 



^*We would see this non-locality in the transverse direction by transforming to momentum space which 
would place a momentum conserving delta function in the transverse plane between the operators. 

^^Formally, the matching coefficient C is a function of SCET label operators but they only appear in a 
Lorentz invariant combination reducing to m|. 
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The hadronic states fix the sum over cohinear directions to be along the protons' initial states, 
and we have made use of the color singlet constraint on the hadronic matrix elements. Finally, 
we have assumed that the so-called Glauber mode does not contribute to the physical cross 
section. Proofs that these modes don't contribute have been given in the more traditional 
approach to factorization [27, 28] but within the EFT, where it is known that Glaubers may 
contribute at the level of amplitudes [26, 29], a proof is still lacking. 

5.2.3 Factorization 

Given the factorized matrix element we now multipole expand it to generate an expression 
for the cross-section which scales homogeneously in the power counting parameter 

_ f d^Ph ,^..+ ^2 _2n/, 1. P^\;,/„2 |2^ 



dp\dy 8v^S{N^-l)J (27r)4^ ' '^'^ V 2 p- 

A{2Trf j d^xe-'^-P>^H{mh)f'l'g,p{Q,x+,x^)f^gip^,{x-,Q,x^)S{Q,Q,x^) (5.19) 

which is valid at leading order in A. We have defined the functions, with spin averaging 
implicit^*^ : 

5(0,0,fx) = (2^)2(^2 _i) (0|'^r(^)^f(^)'gn^(0)'g^"(0)|0)> 

We Fourier transform now to express the factorization theorem directly in transverse momen- 
tum space: 



f±g/pi^^P±) = • Pn){Pn\[B^^{0)5{pnZ - Vn)6^^\p± " A)B^m]\Pn) , (5.22) 

S{0,0,P±) = (jv2^i) (0|^r(0)^f(0)^^(P± - V^)S',%0)S^\0)\0) , (5.23) 

where V is the SCET label-momentum operator. Using the on-shell constraint for the hogs 
and the rapidity delta function, we may set p^ = nihe^y, then, in terms of the momentum 
space TMDPDF^i, we get: 

da TrC?H(mh) f ,0^ /" ,2^ /" ,2^ t-/ 2 1^ -* ^ i2n 

d^y = 2v^S^{N^ - 1) Jdpi±Jdp2^Jd ps^6{p^ - \p,^ +m + Ps±\ ) 

/r./p(^«-^^^i^)/^./P..(3«^^^2x)5(p.J. (5.24) 



^"in what follows, we will denote both the function and its Fourier transform by the same symbol. 
^'^TMDPDF with analogous definitions has been discussed intensely in various contents [55-62]. 
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5.3 Renormalization of Higgs pt Distribution 

Note that to renormalize the soft function and the TMDPDF in transverse momentum space, 
one must adopt a scheme hke 't Hooft-Veltmann or CDR2 [63]. This ensures that the bare 
operator has an integral number of mass dimensions. In 't Hooft-Veltmann, any observed 
degree of freedom is in four space-time dimensions. Any loop momenta, spin averages or 
sum, or internal polarization sums in loops are performed in d space-time dimensions. This 
implies for the TMDPDF f^g^p{z,p±), the polarizations fi and u are in four space-time 
dimensions, and the transverse momentum p± is in 2-dimensions. The proton is spin averaged, 
and thus this is performed in d-dimensions. In CDR2, one allows the polarizations fj, and i' 
to be continued to d dimensions, while keeping the observed transverse momentum in two 
dimensions. 

5.3.1 Renormalization of the TMDPDF 

In perturbation theory, the bare TMDPDF suffers from infra-red, ultra-violet, and rapidity 
divergences. We can renormalize the ultra-violet and rapidity divergences, while the infra- 
red divergence is part of the matrix element. Non-perturbatively, this infra-red divergence 
is cutoff in the hadronic matrix element. For perturbative values of transverse momentum, 
which we will focus on in this paper, we can match the TMDPDF onto traditional PDF's 
and other higher twist hadronic matrix elements (c.f., (5.40)). In the matching procedure, 
the infra-red divergence is canceled, leaving a finite matching coefficient. 
The relation between the bare and renormalized TMDPDF is given as: 

/±5/p(^'^-J-) = Z^^{n,u}/u,p±) (g)± f^llJp{z,p±,fi,uj/u), (5.25) 

where u is the large momentum component carried by the struck parton, the superscripts B 
and R mean bare and renormalized respectively, and we make use of the notation: 

9 ^± f{p) = I ^9{P± - Q±)fiq±)- (5-26) 
In this space we normalize the identity operator as follows 

I^{27r)'6^'\k) = I ^Zj^{k-k')ZfJk'). (5.27) 

The anomalous dimensions of the TMDPDF are then given by 

7^ (P±, /«) = -{Z^^y'0u^Zf^ {p^,fi, uj/u) , 
l-fl^{p±,l^,io/u) = -{Zf^)-^iij-^Zf^{p^,^i,Lo/v). (5.28) 

Notice that 7^ must necessarily be proportional to 6^'^\p±) since hard anomalous dimen- 
sion must be diagonal in pj_ space and the sum of the anomalous dimensions must vanish. In 
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principle there could be plus function dependence on in the TMDPDF //-anomalous di- 
mension which could cancel with the soft function contribution, but given that the TMDPDF 
and the soft function both are renormalized at the same /j, scale, were such contributions to 
the anomalous dimensions to appear they would not contribute to any running. Thus from 
here on we will drop the implied p± dependence in 7^ for both the TMDPDF and the soft 
function. 

The renormalized function then satisfies the RG and RRG equations: 

IJ'j^f±gJp{z,p±,n,uj/iy) =-/l^{iJ,,Lo/i^)f^^Jp{z,p±,n,uj/iy). (5.29) 

5.3.2 Renormalization of the Soft and Hard Functions 

The treatment of the bare soft and hard functions follows in the same way, 

H^irrih) = Z^'ifi, ml)H''{mh, /x) , (5.30) 
SP{p±) = Z^in, ixjv) ®^ S\p^, /X, iijv) . (5.31) 

The anomalous dimensions of the soft function are 

7u{P±,fJ') = -{Z^)~^ ^± y^Z^{p^,iJi,oj/u) , 

7^(/x,M/z^)l5 = -{Z'^r^ ®± ^l^^Z^{p^,^i,UJ/v). (5.32) 

The renormalized soft function then satisfies the RG and RRG equations: 

J^^5^(pj_, 11, 11/ v) = J^{p±,l^) 0± S^{p±,ii, ji/v) , (5.33) 

f^-^S^ip±, fi, nlv) = 7^ (/X, u/u)S^{p±, /X, fi/u) , (5.34) 

where again the /x running can only change the large momentum component. The hard 
function has anomalous dimension: 

if^iruh, /x) = -{Z")-^fi^Z''{fi, mh) , (5.35) 

and satisfies the RG equation: 

f^-^H^'imh, /x) = 7^ (m,,, li)H'^{mh, /x) . (5.36) 

As in the case of the Sudakov form factor we have a set of constraints which the anomalous 
dimensions must obey. The independence of the physical cross section from /x and v gives 

= l^{mh, /x) + 7^^(iv, /x) + 27/^ (w/i/, /x) (5.37) 
= 7;^(/^,P±) + 27^(m,P±)- (5.38) 
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This provides an important consistency check on the calculations of each sector. Furthermore, 
we also have the commutativity of the /x and v running, leading to 

Here we have used the linearity of /i-anomalous dimensions in its logarithmic term and its 
relationship to the cusp anomalous dimension. 

5.4 TMDPDF 

In calculating the transverse momentum dependent PDF, it is useful to consider its matching 
onto the PDF. This will allow us to separate the ultra-violet from infra-red divergences. So as 
long as p_L > Kqcd^ we can perform this matching so that the non-perturbatively effects lie 
in the PDF and its power corrections. The matching onto the PDF is similar to the matching 
of the so-called beam function in [64] 



R I dz' ( 



+ (^ + f^)^X2.A.(^A',pl)}/f/p(^') + 0(^), (5.40) 

where the sum is on species of partons, and the gluon PDF^^ is defined as 

fg/p{z) = -zn-pne{z)g±^APn I [B'J'^{0)6{zn-pn-P)B^\{0)] \ (5.41) 

We adopt the mostly minus metric such that conventions that p^Pj_g±ai3 = —p]_- We make use 
of the 't Hooft-Veltmann scheme, so the external transverse momenta remains in 2 dimensions, 
as do the external polarizations on the operator (the free Lorentz induces). The scheme choice 
is advantageous, as it allows one to renormalize the operator directly in p± space. At tree 
level in perturbation theory we have for the TMDPDF and its matching coefficient to the 
PDF: 

^it/,(^'P^) = '^(l-^M^'Hp±) (5.43) 
lfl/,(^^P±) = (5-44) 
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Figure 7. Diagrams which contribute to the TMDPDF at one loop. Diagram (f) corresponds to 
quark mixing that does not contribute to the one loop running. Diagram (e) vanishes in Feynman 
gauge. 

5.4.1 One-Loop Calculation 

At one loop we have for the sum of diagrams in fig. 7(a)- (e), 



J±g/g ^^^P^) 



2e 



9t 



2 J {2'kY 
2 



5^+\k^) 
{k+Y 



^^^^'^ 2{k-Y+r, j+il ^\uj_-k-) 



(5.45) 



Note that pj_ is strictly two-dimensional, and /c_l is in 2 — 2e dimensions. The momentum 
conservation delta function constrains only the components of k±_ that overlap physical space. 
a;_ is the large component of the light-cone momenta of the incoming parton. Performing the 
integrals we get the bare TMDPDF at one-loop: 



g^CATi-' r(l + e) 



2e 



(2^) 



3-2e 



2 (p^)(i+^) 



-2z + ^ + 



u 



{l-z) 



z"^ p'j_ 



(5.46) 



The virtual diagrams vanished in dimensional regularization. Their effects however are non- 
vanishing and will be accounted for when we match onto the PDF below. Plus function 



Quark mixing is irrelevant for the purposes of this paper. 
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expanding in rj and inserting the MS factor gives: 



J±g/9 l^'P^J 



(27r) 



3-2e 



2 



2 
(1 



5(1 



7/ 

Z) 



2e' 



(1 



0/3 
5± 



(1 



/ p1p1_ g±_ 



Where we have written the expression in terms of the gluon sphtting function: 



Pgg*{z) 



2z 



+ 20(1 



rl 



+ z{l-z) 



(5.47) 



(5.48) 



5.4.2 Matching and Renormalization 

As we used dimensional regularization to regulate both the UV and IR of the TMDPDF, we 



match onto the PDF which allows us to distinguish the — and — poles. Since the effective 



^ and ^ 

<^UV <^IR 

theory containing the PDFs (where the scale set by the transverse momentum has been 
integrated out) must have the same infra-red divergence as the effective theory containing 
the TMDPDFs, any IR poles will cancel in the matching procedure. Any poles left over must 
be UV in origin, and are removed by the Z-factor in the high scale theory. The effects of 
the virtuals, which vanish in dimensional regularization, will arise via the conversion of an 
IR pole into a UV pole. To determine the matching and Z-factor to one loop, we derive a 
relation by expanding equations (5.40) and (5.25) to order a^: 



Z 



(1) 
/x 



c:(-px)+tf„r''(-,px, 



-g/g 



al3 



We split the bare TMDPDF into divergent and finite pieces, / 

idiv; 
-g/g 



(IB) 



^(Idiv) _^^(lfin)^ 



«^^^.r(l + e)-'^'"''^'' '^^^ 



27r2 



(p2)(l+e) J] 



Z) a/3 

-9a: 



asCA S^'^\p±) a/3 f Pgg*{z) 
OsCa 



6{1 - z) In 



(5.50) 



a(5 

P±\9± 



In 



UsCa 1 



vr 



n J \ Pi 



(5.51) 



where we have written the expression in terms of plus distribution £^ 



li^ In" ' 



whose definition and properties are collected in the appendix (F). Given that the renormalized 
PDF at one-loop only contains IR divergences. 



1 

e 27r 



1 



CAPgg*{z) + -PoS{l - Z) 



(5.52) 
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we compare e divergences in eqn. (5.49) to obtain the renormalization constant Z^^ at one 
loop 

(5.53) 

Comparing finite pieces in (5.49), we get the matching coefficients 



= -2"^^£o(/x,Px) . (5.54) 

Note that there is no real singularity as p_L — )■ in I_l2 since the traceless tensor itself 
vanishes in that limit and hence the plus prescription in i^^^gjg be dropped. There are 
finite contributions to the matching from the quark PDF's. These contributions do not effect 
the running of the TMDPDF and hence for the purposes of this paper, we neglect these 
contributions. 

The one loop anomalous dimensions can be calculated from ^j-^"*, 

(5.55) 

ll^ ('^' ^) ^ -SttOsCa^Co (^/x, . (5.56) 

Notice that as long as the scale is taken at the low scale this anomalous dimensions captures 
all the physics at NLO. That is, we need not integrate the [i anomalous dimensions (7^) as 
in (4.30) to calculate since the difference involves no large logarithms. Using the identity 

we see that the results for the anomalous dimensions are consistent with RRG commutativity 
(5.39). 

5.5 The -pT Dependent Soft Function 

The tree level soft function is simply: 

5(0)(p^) = 5(2)(^^)_ (5.58) 

The calculation of the one-loop soft function proceeds as: 
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Using the delta functions and performing the transverse momentum integrals gives 



2(gAg2(7r)-^r(l + e + r?/2) z^V^ 



dk- 



-n 

k- 



(5.60) 



The last integral contains the rapidity divergences coming from both the large k and k~^ = 
limits. This is critical for the soft function in order to cancel the rapidity divergences 
found in both of the jet sectors, as well as consistent with the isotropic nature of soft radiation. 
Finally for the bare soft function we have (including the MS factor): 

2CW e^-'-u^fi'^ r(i + 6 + 2) - 2)r(2) 



(27r)3 (p-2M- 



-I r(i + i) 



vr 



(5.61) 



Since the soft function is IR safe, expanding the p_L-space soft function, we have for the 
divergences and renormalized part: 



2tt;^r(l + e 



1 e''^-B^2e 



1 In^ 



.P±) 



vr 



24 



^5^^\p^) + lnC^)Co(f.A)-C^(,, 



fi J 



P± 



(5.62) 



where w is the book-keeping parameter that tracks the number of eikonal vertices, see section 



(4.1). The plus distribution Cr 



4 In" f C 



Then the anomalous dimensions at one loop order are: 

7„ = —2 m 



is defined in appendix (F). 



vr 



p± 



(5.63) 



This verifies the constraint (5.39) for the soft anomalous dimensions, and on comparing with 
(5.55) we see that we correctly reproduce the constraints 



7^ + 27^^=0. 
Using the result for the hard anomalous dimension. 



(5.64) 



in conjunction with (5.55) (and its partner with lo. 
constraint 



7^+7)^ + 27/^ 







(5.65) 

UJ+) and (5.63) we find that the 

(5.66) 



is also satisfied after making the identification uj^uj^ = Q^. 



- 32 - 



5.6 Resummation of Rapidity Logarithms and Handling Undesired Singularities 



To calculate the resummed cross-section we must evolve the soft function in v up to the jet 
scale z/j ~ as shown in fig. 4. For this purpose we need to solve the Z/-RGE and obtain Vs at 
NLL. First we solve in 6-space, and transform back to present the p_L space solution. We will 
show that a naive solution will lead to a well-known unwanted singularity [65] of the classic 
CSS result [44] . Then we will discuss a careful solution to avoid the undesired singularity. 
Solving (5.33) in impact-parameter space gives: 



5(6, v) = Exp 



vr V 4 / Vz/Q 



where 



In 



vr \VQ, 

Performing the inverse transform gives the resummed soft function in p_L space, 



(2vr)^ 



_L -p_L,M,l^0j , 



where, 



Ks(#±,/",z^, t'o) = 27ry fi6 6 Jo(6|pj_|)^^ ^ 

-20^,7^ ^(1 -tJ ,) 1 



47re' 



Now, the NLL cross-section for p_L > Kqcd is given by 



(5.67) 
(5.68) 

(5.69) 

(5.70) 
(5.71) 



da 



dp]_dy 



Pi >o 



2v^S^{N^ - 1 



-UH{ml,fx^ 



S \e-2^=T^=r(l 



2ml 



1—U)a 



(5.72) 



Note the singularity at Ws = 1 in eqn. (5.72) and (5.71). This is unavoidable because cj^ > 
and is typically ~ 1. Thus relation (5.72) is not useful for phenomenology, but is useful for 
generating the fixed order logs at higher orders in perturbation theory. 

This singularity arises due to the naive inverse transform performed in eqn. (5.70). Note 
that integral in (5.70) gets a singular contribution from h <^ l/p± ~ 1//^- In particular 
when uJs = 1, the integrand goes like 1/6 for small 6 and integral diverges. This shows up as 
singularity at cu^ = 1 in (5.71). This is a UV problem since 6 is small, and is an unexpected 
situation because only impact parameters of order l/p± are expected to contribute to the 
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inverse Fourier transform. Therefore, care must be taken at this step to avoid contributions 
from the region b <^ 

By making a choice for the scale i^o = 1/b in Us, before performing the inverse transform 
in eqn. (5.70), the b <^ l/p± region is exponentially suppressed, removing the singularity^'^. 
One would typically choose vq ~ p± in the p± space NLL cross-section, to ensure all the large 
logarithms are resummed. Making a choice vq = 1/6 is an equivalent one up to higher order 
effects in resummation^^. With this choice we have 

Vs{p±,f,,u) = 27rJ^ dbbJo{b\p^\)[^^ ) , (5.73) 

and soft function is given by 

S{p±,l^,y) = j -^;^Vs{p'j_,n,u)S{p±-p'j_,^i) , (5.74) 

where 

PCX) 

S{p±,fi) = 2tt / dbbJoib\p±\)S{b,i^,uo = 1/6). (5.75) 



^0 

It is difficult to obtain a closed form expression of (5.73), but it is certainly implementable 
numerically and is free from the undesired singularities. We leave this implementation for a 
future work. Since the scale fo does not appear in the coupling, the choice of scale setting 
does not involve any Landau pole, as /U is left arbitrary in the transforms. 

A similar approach was taken in [40] to cure the singularities, however, there the singu- 
larity was interpreted to be some indication of non-perturbative physics. In contrast, it was 
argued in [65], the problem is completely perturbative: it is solvable by a rearrangement of 
sub-leading terms of the resummed series, and occurs in regions where as is still perturbative. 
We agree with this and have verified that the problem simply arises from a naive inverse trans- 
form that includes contributions from a UV region inappropriate to the soft matrix elements 
being resummed. 

5.7 Fixed Order Cross-section 

Substituting the tree level matching onto PDFs, and integrating over p± and y we obtain 
properly normalized integrated cross-section at LO, 



2v^S^N^ -I) 2ml 
This agrees the leading order cross-section given in [66]. 



dzidz26{m\ - ziZ2S^ fg/p{zi)fg/p{z2). (5.76) 



■^^We thank Wouter Waalewijn for making this suggestion. 

^^Typically we can only argue ~ 1/b but the arbitrariness is still captured by varying v ~ in a 
reasonable range, thus fo — l/b is justified. 
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For the next to leading order transverse momentum spectrum we get 




where g ® h{z) = ^g{x)h{z/x). Note that an overall factor of vr comes from the angular 
integrations in the transverse momentum convolution variables. This agrees with the result 
found in [67]. 

5.8 Comparison to Previous Resummation Formalisms 

Working within an EFT formalism, a result nearly identical to (5.72) was derived in [40] at 
NLL. However, the rapidity logs were exponentiated by hand after summing over sectors^^, 
without introducing a new scale. Thus it is not clear how the scale dependence arising at 
higher orders can be tracked. As explained in detail in the conclusions of sections (4.3) 
and (6), there is a fundamental ambiguity in the exponentiation of the rapidity logarithms, 
since there is freedom in choosing what goes into the exponent one is free to include sub- 
leading logarithms in the resummation power counting in the exponent or the low-scale matrix 
elements. In our formalism this ambiguity corresponds to the choice of matching scale. 
Varying this scale shuffles sub-leading contributions into or out of the matrix element. Hence 
the residual u dependence (which is not exponentiated) of the resummed cross-section, can be 
included in the theoretical error of our prediction by varying v in the same way one varies fi to 
get a handle on errors form sub-leading term in traditional RG calculations. Such an analysis, 
at least in the context of Higgs tranverse momentum distributions, is presently absent from 
the literature 

It is also worth noting that our calculation of the transverse momentum distribution 
is distinct from the work [40]. Therein, they do not include the soft mode of the effective 
theory. It is stated that the "soft mode" cancels because the typical soft momenta is order 
(A^,A^,A^). We would call this mode ultra-soft, which is the relevant mode for SCETj. The 
use of the analytic regulator renders integrals in the actual soft function (having modes 
Q{X, X, X)) scaleless and hence zero. While technically correct, this method seems to obscure 
the physics, given that soft radiation clearly plays a role at small p±, and thus is must be 

■^^This sum over sectors is necessary due to the regulator chosen to give meaning to their low-scale matrix 
elements. Given the regulator implemented in [40], each sector is not well-defined and nor renomalizable, but 
only the combination is. 

■^^Such an analysis can also be performed in the CSS formalism using Collins most recent definition of the 
TMDPDF. For a discussion see [68]. 

^'^The traditional analytic regulator breaks the eikonal identities needed for exponentiation of the soft func- 
tion. This complicates the claim made for factorization, though this breaking of eikonal identities has been 
remedied in a more recent the paper [69]. 
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hidden in other sectors. Finally, the implementation of the analytic regulator in [40] renders 
the expanded results in the two coUinear sectors different in form, even though the operators 
in the effective theory look identical, and thus one can not define universal TMDPDFs. 

The classic CSS formalism utilized the so-called Collins-Soper equation to accomplish 
the rapidity resummation. In particular, it played a key role in establishing the formal /i- 
independence of the double logarithmic terms in the resummation exponent [7]. Using this 
formalism, a resummed formula for transverse momentum distributions was derived [44]. 
Nonetheless, the classic CSS formalism suffers from a number of mild deficiencies according 
to one of the authors[70] . In particular, the Collins-Soper equation itself suffered from hard to 
control power corrections, and the hard matching coefficient is ambiguous. These issues stem 
from the way the rapidity divergences were regulated in the original Collins-Soper approach, 
where a non- light like axial gauge regulated the divergences. The regulating parameter did 
not cancel between the low-scale matrix elements, so the choice of the axial gauge vectors 
that defined the collinear matrix element also affected the hard matching. Nor was the 
regulating parameter divergences ever removed from the matrix element, so the regulating 
parameter could not be set to zero. This malady remained true in more modern versions of 
the Collins-Soper approach adopted in [19, 71] where the Wilson lines are deformed off the 
light-cone. 

The authors of [48, 49, 72] who have performed the highest order resummations of the 
transverse momentum spectrum to date have used the resummation formula given in [44] . The 
anomalous dimensions and matching needed for resummation were derived by comparison to 
full QCD calculations in soft and collinear limits, not by direct calculation of the low-scale 
matrix elements in the Collins-Soper formalism. From this procedure, it is not clear that they 
can gauge the residual effects of the rapidity resummation, since they do not make use of the 
Collins-Soper equation directly, nor calculate the resummation from the factorized matrix 
elements. 

Recently Collins has improved further upon the CSS approach in [70], fixing the above 
problems. This new method has many similarities to our approach. In [70], the square root 
of the soft function is included in the collinear sector, and a series of soft-bin subtractions 
(explicitly represented by inverse soft functions) are carried out to cancel the regulator de- 
pendence in the sector. This introduces explicitly a In — in the TMDPDF, where p~ is the 
large light-cone momentum. The soft factors remove the rapidity regulator dependence and 
introduce a fi dependence in its place. Thus if p~ is large, a single fi scale cannot eliminate 
all large logarithms. However, through use of the Collins-Soper equation, the extra large log- 
arithms can be exponentiated by evolving the TMDPDF from low energies, where p~ ~ 
to the relevant high energies. In the low energy region, all logs can by minimized by a single 
choice of /x. All logarithms remain minimized when evolving the low-energy TMDPDF via the 
Collins-Soper equation to the high scale where the experiment takes place. This accomplishes 
the same effect as our rapidity RG, but we do not need the low-energy TMDPDF as an initial 
condition (just the PDF). We can take the TMDPDF at fixed energy and transverse momen- 
tum to be renormalized at any |U, z/ point, and evolve to any other point in order to minimize 
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logarithms. Finally, we must mention that there is one important way in which the claims^^ 
[70] are stronger then ours. In particular the proof in [73] allows for the possible contribution 
of Glauber gluons. In SCET an understanding of Glauber gluons in hard scattering processes 
[26, 74] is still lacking. 

Another attempt at defining a TMDPDF within the SCET context was presented in 
[75]. Using a <5-style regulator, and including a square-root of the inverse soft function in the 
collinear sector, they were able to eliminate the rapidity divergences within the bare matrix 
element as in Collins' approach. They do not discuss resummation, nor indicate how, after 
eliminating the rapidity divergence with no auxiliary parameter, they will accomplish the 
specific resummation of the rapidity logarithms. As it stands no choice of /U will minimize all 
logarithms in their TMDPDF if the light-cone momentum is large. 

Finally, it is worth noting that the approaches taken in [75] assumes the equivalence of 
soft-bin subtractions and the inverse soft function. While this is true in many situations for 
many regulators, this is complicated in the case of a soft function depending on multiple 
parameters. This is illustrated below, when we factorize the cross-section for left and right 
broadening. Here the soft function depends on the broadening in both hemispheres, while 
the jet function only depends on the broadening of a single hemisphere. Thus the soft bin 
of the jet sector cannot be the the same as the inverse soft function, as the soft-bin will 
continue to only depend on the broadening of a single hemisphere. Including the square- 
root of the (inverse) soft function to make the jet function free of rapidity divergences would 
introduce dependence on the dynamics of the other hemisphere. The fact that the soft function 
depends on both hemispheres complicates a straightforward definition of a square-root of the 
soft function, especially since non-global logarithms seem to be a generic feature of multi- 
region soft functions [76-78]. Thus it is hard to see how one would generalize Collins' recent 
approach to the TMDPDF to cases such as jet broadening where one cannot eliminate the 
soft function by splitting it up the between the jet functions without inducing dependence on 
the broadening of both hemispheres in each jet function. 

Finally, for other works within SCET transverse momentum factorization [41-43] only 
the hard logarithms from running the currents are resummed. The rapidity logarithms are 
left unsummed. 

6 Jet Broadening 

Event shapes have played an important role in precision measurements of the strong coupling 
as[79]. A generalized event shape for event e~e^ ^ X at center of mass energy ^/s, can be 
defined [80] in terms of a parameter a via 
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We use this terms simply because we have not been able to reconstruct the proof ourselves. 
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where pj_L is the transverse momentum with respect to the thrust axis t of the event, and rji 
is the rapidity of the i'th particle. The thrust axis t is defined by maximizing thrust T [81], 



T close to 1 corresponds to the special case a = 0, e{0) ^ 1 — T and is also loosely called 
"thrust". Another interesting event shapes is the limit a = 1 corresponding to "total jet 
broadening" B [3], e(l) = 2B. The limit e(a) <C 1 isolates events composed of back to back 
jets. In the case of thrust, jets are composed of collinear radiation, and the recoil due to soft 
(ultra-soft in this case) radiation does not affect the jet axis. For jet broadening all radiation 
with parametrically similar transverse momentum can contribute, so the soft radiation of 
order Q(A,A, A) recoils the jet off the thrust axis, where A ~ e(l). In both of these cases 
fixed order perturbation theory will fail when e(a) is small. However, as long as eQ ^ Aqcd, 
we expect non-perturbative effects to be suppressed, though large logarithms of e need to be 
resummed. 

The pioneering work on jet broadening resummations [3] utilized the coherent branching 
formalism [82]. It was later stated [2] that the results in [3] neglected terms due to recoil of 
soft gluons. In this section we will provide a factorization theorem for jet broadening. The 
factorization proofs for angularity observables (6.1) in [83] are known to fail as a approaches 1, 
since there are growing power corrections in this limit. The reason for the apparent breakdown 
of factorization is the fact that in this limit the soft radiation has the same invariant mass 
as collinear radiation and one must change the power counting accordingly to factorize in a 
consistent fashion. Which is to say that thrust can be analyzed using SCETj whereas jet 
broadening necessitates the use of SCETn. 

6.1 Factorization Theorem 

We start with the expression for differential cross section in QCD for broadening e (strictly 
speaking, angularity for a = 1), 



? = ^ E / d'xLl,Me'^-^0\jl''{x) {S{e - e)j^m |0) 
i=Ay'' 

= ^ E i;..('?)E(2^)''^^'^(^-^^)Wl^0)5(e-e)|X)(X|jT(0)|0), (6.3) 



where e is the jet broadening operator that first maximizes thrust for a given state \X) to 
determine the thrust axis i and then measures broadening via, e\X) = '^^iex ^^^l"'^)' where 
p±i is transverse momentum measured w.r.t. t. Here Q is center of mass energy and QCD 
current is given by 




(6.2) 



i=A,V 



X 



J. 



/f (x) = q{x)T'^q{x) 



(6.4) 
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with = 7^ or 757^. L^;/ is the leptonic tensor given by 



L 



L 













levr^aem 






3Q2 





:i-m|/Q2)2+r|/m| 



(6.5) 



■"^ ' (l-m|/Q2)2 + r|/m|_ 

where subscript q denotes the (anti)quark flavor, Qg is the quark charge in units of |e|, Vg^e 
and tig g the vector and axial couplings of the (anti)quark q and the electron to the Z as 
e.g. in eq.(A3) of ref. [84]. Here niz and Tz denote the mass and the width of the Z boson. 

In what follows, by requiring e ~ A ^ 1, we will prove a factorization theorem of the 
form 

^ = HJn ®Jn®S, (6.6) 
ae 

where each function is a vacuum matrix element of operators that depend on either collinear 
or soft modes that do not interact. 

We begin by first matching the QCD currents onto the SCETn currents2^ as follows, 



(6.7) 



where Cn^n2 = C{ni-V,n2-V), is only a function of large label momentum operators. Xni,n2 
are SCET collinear fields while the Sn are soft Wilson lines extending to infinity along the n 
direction. Inserting eqn. (6.7) into eqn. (6.3) we have 

da 1 
de 



ni ,n2 hi ,h2 



(6.8) 



X 



mrnSl rf Xn2 S{e - e)\X){X\xn2Si T'^ Sh, XnJO) , 

where it is implicitly understood that all fields are evaluated at space-time coordinate x = 0. 
Now we decompose the complete set of states as = J2h X^x^^ \^n){Xn\ where Xn 

are states with thrust axis along h. Here, the states that may have an ambiguity in the choice 
of the thrust axis, without loss of generality, can be associated with either h. For q^ = {Q, 0) 
we have 
da _ 
de 



'^l,"2 ni,n2 h 



n2^h\h2 



(6.9) 



2Q2 



-^^1/ ^ni^2^fiih2 (0| XniSli-^ Sn2 Xn2 Xn2'5'^2 "^"i 1^-* 



In principle we should include ultra-soft Wilson lines Yn's as well, but since ultrasoft modes do not 
contribute to the process the ¥„ Wilson lines will cancel. Therefore we drop them already and resort to the 
phrase "SCETn current". 
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where we identify the term in braces as the broadening projector X^, which simphfies as 
Xn = Y,^iE-Q) 5(3) (P) 6{e-en) \Xn){Xn\ 



= ^6{E-Q) 5(3) (P) S{e-en)S^,oS^^,o \Xn){Xn\ 
= Y,^iE- Q) 5(3) (p ) 5(e - in) 5^,0 5|^,o I^X^I 

X 

= 5Qn-P+^n.P-Q) ^Qn-P-^n-p) 5(2)(P^ ) 5(e - e^) 5i^,o 5|^,o , 

= 2 5(n-P-Q) 5(n-P-Q)5p^o S{e - in) 6^^\¥n^) 5^^\f^) {^j <fkr^ 

= 2 5(n-P-Q) 5(n-P-Q) 5(e-eft)5(2)(P^)5(2)(P^) {^j (fK^ , (6.10) 

where and are defined as 

K^\X) = Y,0in- Pj)Pj±\X) 

fn±\X) = Y,Gi-^-Pj)Pj±\X)- (6.11) 

In the second Hne we have simply used ^i^^ol-'^n) = \Xn) and S^^^lXn) = \Xn) owing 
to the property of the thrust axis that total transverse momentum in each hemisphere de- 
fined by the plane perpendicular to the thrust axis is zero. This form of the broadening 
projector implements the kinematic constraints imposed by the choice of thrust axis, since 
zero-transverse momentum flow in each hemisphere defined by h along with small broadening 
w.r.t. fi uniquely fixes the thrust axis to be n. Hence in the third line we have promoted the 
state \Xn) to a generic state \X). In the fourth line of eqn. (6.10) we have summed over the 
complete set of states. For convenience in factorization with continuous labels, we turn the 
Kronecker-5s to continuous Dirac-5s, with the general relation: 

5p, o = '5^'^(Pi) / d'^r-±, (6.12) 



where / d?kj± is simply the area of the label-transverse momentum unit cell^*^; refers to 
the residual momentum. In the final step we have used q = ^XT±+%±fi ~ ^^fi ~ ^' 
e <C 1, constraints put in by all the 5- functions ensure that the broadening projector selects 
the dijet states with the thrust axis h. 

Given the properties of the broadening projector for small broadening, the light-cone 
directions of the collinear fields must be within a small cone (of radius A) about the h- 
collinear directions. Fields whose directions are an order 1 or more displaced from n must 



^°The Kronecker-(5 assures that the discrete-label momentum is zero while the Dirac-(5 ensures that 
continuous-label momentum is zero up to order A^, i.e. all momenta belonging to the unit cell are considered 
to be zero, hence the area of the unit cell appears as the proportionality constant in the identity (6.12). 
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produce too great a broadening, and so are excluded. Therefore by making a parametrization 
transformation [5], we can set the directions of the colhnear fields to be exactly parallel or 
anti-parallel to n. With further constraints from quark number conservation we only have two 
choices: the quark jet is either along the n or —h direction. Since the observable is symmetric 
under charge conjugation, the two are equivalent, and we get 

^ = ^^M- 2 CnnCln (0| XnSi Sn Xn XnSl T\ Sn Xn |0) (6.13) 



2 \Cnn\^ (0| XnSi Tf Sn Xn 5{n ■ P - Q) 6{n-P- Q) 



lb 

x5{e-en)5{K±)5{fr±) XuS^T^i SnXn |0) ( / <fk 



r± 

where we have used eqn. (6.10). The matrix element in the last equation does not depend 
upon the direction h but only on the large labels Q and broadening e. The choice of n inside 
the matrix element is only representative of a light cone vector necessary for calculation, but 
these calculations would yield identical results for different n. Also, Cnn is only a function of 
boost invariant n-Pn-P = s = Q^, so we define the hard function independent of the light 
cone direction, H{Q'^,ij,) = |Cnnp. Therefore we can safely factor out the matrix element out 
of the sum over h and can write 

% = ^^^"^T^ H{Q',fi) (0| XnSi rf Sn Xn 5{n ■ P - Q) Sin-P- Q) (6.14) 



x6ie-en)6iFn±)6iFnL) XnSiV^SnXn |0) hV / d\r± 




where h from here on is a fixed vector, say z, i.e. = (1,0, 0, 1) and = (1, 0, 0, —1). Now 
we use [85] 

/ d^k± X] = 2^ / d'^kr± = ^ f dn = TrQ\ (6.15) 

cones n 

where k± is the label-momentum and kr± is the residual momentum. The cones subtend an 
angular area of order 1, while h directions label cones of order A. Then we achieve: 

^ = i27r)\HiQ^^)Ll,iq) (6.16) 

{0\XnSl rf Sn Xn 5{n ■ P - Q) 5{n-P- Q) 5{e - en) 5^''\FnA_) 5^^\fnL)XnSi rr 5^ xn|0>. 

We have two choices for lepton tensor which correspond to = 7^ or 757^ in the hadron 
matrix element. We can simplify the lepton tensor by noting that Xn {<i-> lb4)Xn = for 
g_L = 0, thus we can safely replace P^^^iq) with U{Q'^) g^y. We now have 

^ = (27r)V/7(g2,^)LXQ') (6.17) 
(0| Xn4rf 5„Xn5(n-P-g) 6{n-P-Q) 5{e - e^) 5^^\¥nA_) 5^^\¥nA_)XnSiV^SnXn |0). 
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We are almost ready to factorize. The operators e^, and still mix n-collinear, n- 
collinear and soft sectors. We take care of this by considering following identity operators 

1 = J d'k^ d'k2L <fk'^ <fk'^ <5(2)(pW _ fc^) 5(2) (^^5) _ 5(2)(pW _ 5(2)(^(s) _ 

(6.18) 

and 

1 = j den den des 5{en - ef ) 5{en - ef) 5{es - 4'^) , (6-19) 

where operators with superscript (c) has the same action as their parent operator on n- 
collinear particles and fields but they give zero for all else. Similarly (c) for n-collinear and 
(s) for soft. We now insert (6.18) and (6.19) in the expression (6.17). Since there are only 
n-collinear, n-collinear and soft sectors contributing to (6.16) we must have wj:^ = ^nL, 

+ = I^lL and Cn = e^'^'* -|- e^^^ -|- e^f^ . Using this and integrating over k'^ and A;^ , we 
have 



{2t,)\ H{Q^,^j)V{Q^) j den den des J d^k^d^k2L (0| Xn4rf 5„Xn 



da 
de 

x6{n-P-Q) 6{n-P- Q) 5{e - en - en - e.)(5(e„ - ef)6{en - ef )<5(e, - 4'^) 

y<5^^\¥;^-kn_)5^'){f;^-k^)5^'\F;^ + kn_)5^'){fi^+k2^) XnSlr.SnXn |0) . 

(6.20) 

The hard work is done, now we can simply factorize because each operator acts only on either 
sector (note that n • P and n • P get leading contribution only from n-collinear and n-collinear 
sectors respectively; in other words multipole expansion ensures that we can safely assume 
n • P acts only on Xn held and so on) . After using color conservation for the collinear matrix 
elements and the Fierz transformations, 7ai,7Mcd -(|)ad(|)c6 and {l^Y)ab{l^liJL)cd 
-(|)ad(|)cb , we have. 



Nci — jH{Q'^,n) I dendendes 5{e - en- en- es) d'^ki±d'^k2i_ 



de 



^ (0| Xn d{n-P- Q) Sien - en) d^^HP± - k^i) ^ Xn |0) 



(2vr)3 



tr (0| IxnKn-P- Q) 5{en - en) 6^^\P^ - k^) Xn |0) 



^tr(0| SiSn 6^^\Fn^ + k^) <5(2)(]P^ + A;^) 6{es-en) SlSn |0) , (6.21) 

where we have also reduced transverse momentum and broadening operators to their re- 
spective parent operators (with a further reduction of transverse momentum operator in the 
collinear sectors to label momentum operators). Trace is over both color and Dirac indices. In 
eqn. (6.21), we have a factorization theorem in which the third and the fourth lines represent 
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the broadening jet function composed only of collinear fields, and in the last line we have 
a matrix element composed of the soft Wilson lines only. Written compactly, factorization 
theorem is 

da 



aoH{Q'^,fj,) j dendendes 5{e - Cn - en - Cs) j d^k^d^k2i_ 

Jn{Q, en, klj_) JniQ, Bfi, klj_)S{es, ki±, k2±) (6.22) 
where ctq is the Born cross-section given by (see, for example, appendix A of ref. [84]) 



3Q2 



(6.23) 



(l-m|/Q2)2 + r|/m| 

A straightforward generalization of this result gives the factorization theorem for the left and 
right broadening [86] 

da 

deideR 

MQ, e„, kl^) JniQ, en, kjx) 5(ef , ef , kj^,kl^) , (6.24) 



aoH{Q'^,fi) J dendende^ def 6{eR- en - ef)6{eL - en - eg) J dk1i_dkli^ 



where only the soft function changes 

2 

5(ef , ef , kj^, kl^) = ^tr (0| 4 5„ + j,^) ^^^R _ g^) 

5^''\fni. + k2i_) 5{e^s - eft) S], Sn\Q). (6.25) 

Here e^'^ are defined as e^'^|X) = {Y^^^-^9{±.pi'j,)\pix\/Q)\X) . This simply comes about by 
changing the identity insertion (6.19) appropriately. Note that in this form each function only 
depends on the magnitude of the transverse momenta which is particularly convenient, thus 
we have changed the overall integration measure appropriately. We present the bare definition 
of the jet and the soft functions in the next section where we discuss their renormalization. 

6.2 Broadening jet and soft functions: definition and renormalization 

The naive definitions of the jet and the soft function obtained in the previous section contain 
unregulated rapidity and UV divergences. We will regulate the UV divergences in dimensional 
regularization as usual and for the rapidity divergences we will adopt the regulator prescribed 
in this work. For fermion free abelian theory we can put the rapidity regulator into the 
collinear and soft Wilson lines and for the non-abelian case we follow sec. A. 2. Bare quark 
jet function is given by 

(2^)3-2^ . M 



J^--^%Q,en,kl) = ^ tr(0|^Xn(0)<^(n-P-Q) 5{en - en)5^^-^'\V^ - k,_)xn{m) , 

(6.26) 



where V is the standard SCET label operator and here we are working with continuous labels. 
There is analogous equation for bare anti-quark jet function but its functional dependence is 
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the same as the quark jet function so it is simply obtained by replacing e„ — )• Cfj. The bare 
soft function is given by 

5bare(ef , kl^, kl^) = ""^^^^I^Z^tr (0| 4(0) S„(0) 5^^-'^\^^ + k^) 5(ef - ) 

5(2-2e)(p^ + A;^) 5{e^ - 4) St (0) 5^(0) |0) , (6.27) 

where and were defined in eqn. (6.11). For both the jet and soft functions all fields are 
evaluated at x = 0. Note that after accounting for the dimensions of the bare fields, bare jet 
and soft functions are integer dimensional objects, as required for operator renormalization. 

In our formalism rapidity divergences appear as counter terms just like the UV diver- 
gences and therefore can be renormalized away via renormalization constants. Renormalized 
quantities can be calculated as usual 

r'^-{e,P-fi,u/Q) = I de' I 0^Zj\e-e',{k-k'f;f,,u/Q)j''''''{e',k'^Q) 

S^'^''{eR,eL,p',f;f^,i^) = J de'^de'^J '^^'''(4, g'^) (6.28) 

^S^i^^R - ^R^f^L - e'i, {p-pf, {q- qf]H,v) , 

where Zj^^ only contain terms that are divergent in r/ and e. Note that after expanding in rj 
and e all vectors are Euclidean 2-vectors, hence we have a 2-dimensional vector convolution 
only. Henceforth, in this section, we drop the subscript _L on 2-vectors. Zj^s follow the 
standard constraints 

Ijie,k) = (27r)V2)(fc)5(e) = J de' J ^Zj\e-e',ik-k'f)Zjie',k'^) (6.29) 

IsieL,eR,p,q) ^ {27r)^S(-^\p)5(^'\meR)S{eL) = J de'^de'^ j ^-5(6'^, ei,p'^ g") 

Zs^{eR - Cr, ei - e'i, {p-p'f, {q- qf) ■ 

Using the consistency condition and the fact that the bare functions do not depend upon 
the renormalization scales // and v we get the RG equations in /i and v. The renormalization 
group equations for the jet functions are given by 



^ldJ 



ren 



^)i{^Ji,u/Q)J■'''''{e,k^■^^,u/Q) (6.30) 



de' / ^J^{e-e',{k-k'f■,^^)J^^^e',k"■,^^,u/Q), 



du J J (27r) 

where the anomalous dimensions are obtained via 

72 D 



x^^{fi,u/Q) = -l de' I ^Z-\e-e',{k-k'f)^Zj{e',k";^,u/Q) (6.31) 
^^{e,P;f,) = - I de' I ^Zj\e - e',ik - k'f)^Zj{e',k'''; f^,u/Q) . 
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For the soft function, RG equations are similar with more variable dependencies, therefore 
for brevity we only show relevant variables and represent convolutions via (8). We have for 
RG equations 

,. cRena 

^^^^ = 7;!(M,^^^)5-'^(ei^,ei,^r^g^M,z.//i) (6.32) 

and for anomalous dimensions 

Is x-f^^{f,,u/Q) = -Z^^(^^^Zs{...;fi,u/f,) (6.33) 

75(eR>eL,p^^^;/x) = -Z^^ O J^^si- ■ ■ ; fJ-, / fJ-) ■ 

We emphasize that //-RG equations do not involve convolutions and 7^ does not have any 
kinematical dependence. Reasons for this were made clear in section 5.3.1. 

We have a few consistency conditions on the anomalous dimensions. Firstly, sum of the 
/x-anomalous dimensions in the IR sectors (jets and soft) should add up to the negative of the 
hard anomalous dimension, that is 

2 

7;f ill, In ^) + 7^ (/i, In ^) + 7^ (/U, In ^) + -ff! (fi, In ^^) = , (6.34) 

where we have explicitly shown the logarithmic dependence in the anomalous dimensions 
which can only be linear. Note that in c.o.m. frame Q"*" = Q~ = Q. Secondly, we have a 
consistency condition for i/-anomalous dimensions analogous to eqn. (4.28) 

u4- [{ZjZj) ^Zs]=0 



dv 



fhi^ {e,f) + Ij{e,f) 1^ (e', f) + -il{e, e' ,p\t) = . (6.35) 



Lastly, the independence of renormalization scales, /u and v implies that the UV and rapidity 
RGs commute, which then gives the constraints 



_d_^j7_j ^d^ 
dfj, ^ dv 



i^— = = r^sp ij, 

d „ d 



/i^Ti^ =%^^^7m = -2n,spl5. (6.36) 

Here we have used the linearity of /i-anomalous dimensions in its logarithmic term and its 
relationship to the cusp anomalous dimension in 7^. This concludes the formal discussion on 
renormalization and running. We present the calculation and results in rest of this section. 
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6.3 Jet Function Calculation up-to NLO 

Bare quark jet function was defined in eqn. (6.26), which is what we will calculate here. 
Anti-quark jet function is obtained by e/j — >• e/, and — >• g^. For tree level jet function we 
have 



J(0)( 



HeR-\p\/Q), 



(6.37) 



which is all what we need for the NLL cross-section. For NLO-singular cross-section, we only 
need to consider the one-loop jet function at p = 0. Non-zero transverse momentum implies 
presence of at least a soft radiation, which then gives contribution to the cross-section at 
NNLO. We will use ?7-regulator as prescribed in section 4.1 to regulate rapidity divergences. 
To regulate IR and UV divergences we will use dimensional regularization with d = 4 — 2e. 
Virtual diagrams all vanish as they will be scaleless. So we need to calculate only the real 
diagrams and interpret l/x^"*"" as the distribution [;^rT7r]+^^- We use Feynman gauge for this 
calculation. At one- loop, the jet function is simply the sum of all real diagrams^^. We get, 
for p = 0, 



(a) 



(b) 




(d) 





Figure 8. Diagrams contribution to broadening jet function. Diagram (c) is zero in Feynman gauge 
and virtual diagrams (d) and (e) are zero in dim.-reg. 
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(6.38) 



^^this is the plus distribution with boundary at +oo i.e., dx[ _^i^l ]^ = 0. 

■^^We have verified by using gluon mass as an explicit IR regulator that IR divergences cancel between the 
real and virtual diagrams and that we obtain the same results as presented here. 
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Expanding in rj and then in e we get, 
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where we have expressed the finite parts in terms of the standard plus distributions which 
are related to those with the infinity boundary via 
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and obey 
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Since we did not calculate p'^ dependence we only give the renormalized part of the one loop 
jet function which is what we need for the NLO cross-section, 
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(6.42) 



6.4 Soft Function Calculation up-to NLO 

The bare soft function was defined in eqn. (6.27). The tree-level soft-function is given by 

5^ = 6{eR)6{eL)5{p^)d{f) . (6.43) 
At one-loop we regulate rapidity divergences with r] and use dimensional regularization for 




Figure 9. Diagrams contribution to broadening soft function. The virtual diagrams (c) and (d) are 
zero in pure dim.-reg. 
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UV and IR divergences, just as in the previous section^'^. In Feynman gauge only two real 
radiation diagrams contribute to the one loop soft function which are identical to each other 
leaving 



^ n.k n.k 
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Expanding in r/, then in e and combining with the tree- level result we have the bare soft 
function at NLO, 
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As in the jet case, we have written the rapidity divergences in terms of plus-distributions with 
boundary at 00 and the finite corrections in terms of the standard plus-distributions. 
Using eqn. (6.28) we extract the renormalization constant for the soft function. 
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(6.47) 



We have verified that regulating IR divergences with gluon mass at one-loop yields the same results as 
presented here. 
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The anomalous dimensions can be calculated using eqn. (6.33). At one-loop we have, 
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It is a straightforward exercise to check that constraint (6.36) for soft anomalous dimensions 
are satisfied up to order Og- We cannot do a direct check on constraints (6.34) and (6.35) but 
an indirect check comes from calculating the cross-section at order ag with the bare matrix 
elements. We will perform this check in sec. 6.7. 

6.5 NLL Soft Function 

To calculate the NLL cross-section for jet broadening we must evolve the soft function in ly 
up to the jet scale (Q) as shown in fig. 4. For this purpose we need to solve the v-RGE and 
obtain Vs at NLL. Its easiest to solve this equation in conjugate space where we make Fourier 
transform w.r.t. p, q and make Laplace transform w.r.t. ei, en- In the conjugate space 
Z/-RGE reads 
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where bji (bi) is Fourier conjugate to p (q) and tr^l is Laplace conjugate to cr^l- Since 'ji, 
has no explicit dependence on v solution of this RGE is simply 



5(. . . ;/U,i////) = exp( %(. ..;/i) In — ) 5(. . . ; /i, fo//") • 



(6.50) 



Therefore, 



Vs = LF 



-1 



% {TR^TL.bR^bL-.^i) In — 

^0 



(6.51) 



where LF ^ refers to the inverse Laplace and inverse Fourier transform on tr^l and 6_r,l5 
respectively. From Laplace-Fourier transform of the result in eqn. (6.48) we find that 
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which gives 
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where 



uJs = a;NLL(as(/"),z^/z^o) 



(6.54) 



vr vq 

Using eqns. (E.4) and (E.5), we evaluate inverse Laplace-Fourier transforms in eqn. (6.53), 
to obtain the i^-evolution factor Vs at NLL in the physical space, 
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Note that lOs is always positive for us (i^ > fo), so there is no non-integrable singularity in the 
last two fractions and other terms are properly plussed. Convolving this with the tree-level 
soft function gives us the NLL soft function, 



5nll = ilJ-Q) 



2 -2uJs^E 

'2uJs 



r2(a;. 



0{QeR-\p\) 0{QeL - \q\) 

72\l-uJs 



(Q2e2 



R-P 



(6.56) 

This very interesting formula illustrates how the soft function changes when we take into 
account infinite gluon emissions from a quark between the rapidities of order | ln(fo//^)| ^-nd 
I ln(z^//i)|. Resummation of all the large rapidity logarithms require that we choose ~ z^o ~ 
QeL,R and v ^ Q. This would mean that 5nll accounts for recoil from gluon emissions of 
rapidities, \y\ < |ln(eL^ij)|. Note that the unwanted singularities that arose in section 5.6 
do not arise here. Reason being, there is an absolute cutoff on the transverse momentum 
governed by the measured eL,R- Thus no unwanted UV contributions can arise in the inverse 
transform. Alternately, there is no singularity from 6 — )• in eqn. (6.53) for any value of ojg- 

6.6 Jet Broadening Spectrum at NLL 

It is a straight forward exercise to calculate NLL differential cross-section now. NLL formula 
is given by 
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where i^tree = 1 and is the evolution factor for the hard running. It is the same as for 

other angularities (for example thrust) and can be found in, for example App. C.2 of ref. 
[20]. Tree-level jet functions were given in eqn. (6.37) and are simply (5-functions. Taking the 
NLL soft function from eqn. (6.56) and performing straight forward J-function integrals we 
get 
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where we have non-dimensionahzed the last remaining integral and dropped the plus prescrip- 
tion since the integrals are well defined. Performing this integral yields the master formula 
for NLL cross-section 
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where Bz{a, b) = dxx"' ^(1 — x)^ ^, is the incomplete beta function and a;^ was defined in 
eqn. (6.54). Using this result we obtain the NLL cross-section for wide jet broadening Bw 
[3], 
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and total jet broadening [3] 
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6.7 Jet Broadening Spectrum at LO 

We will compute the spectrum at LO using the bare matrix elements to elucidate the can- 
cellation of rapidity divergences and corresponding scale v. This distribution accounting for 
one real radiation is given by 
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which can be computed by using (6.39) and (6.45) along with corresponding tree level results. 
We obtain 
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As expected both the rapidity divergences and z/-dependence cancels out. The sum of the 
UV divergences in the IR sector add up to give the expected form consistent with the hard 
anomalous dimension, i.e they cancel with the UV divergences in the bare hard function 
[38, 87], 
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Using eqn. (6.63) we calculate the singular contribution to the LO total broadening distribu- 
tion 



1 dgLO 
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(6.65) 



which agrees with [3]. 
6.8 Numerics 

In fig. (10) we have plotted the theory cross section and the data [88]. The resummed error 
bands are the geometric mean of the i/- variation and /i- variation. We see that given the large 
experimental error bars the agreement with the data is reasonable. Complete jet function 
calculation up to NLO, although will not change the NLL resummed spectrum, is expected to 
bring down the scale uncertainty significantly. In addition, the NNLL calculation will reduce 
the theory errors further. It is worth noting that by including both u- and ;U- variations, we 
gain a well controlled theoretical uncertainty estimation, while the uncertainty analysis for 
resummation using traditional methods could be ambiguous and may under or over estimate 
the uncertainties as we will discuss in the next section. Here we have not included the theory 
errors due to power corrections. In the small Bt region these are non-perturbative and scale 
as AQCD/iBrQ) and can be expected to be of order 20-30%. 
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Figure 10. Total Jet Broadening at 130 GeV. 



6.9 Comparison to Previous Results 

In our previous work [23] where we introduced the rapidity renormalization group, we pre- 
sented results for NLL total broadening distribution with an unnecessary approximation that 
oversimplified the structure of 7;^. As a consequence, our result presented in [23] was not 
accurate at NLL when power counting the resummation in exponent. The result presented 
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in here in eqn. (6.61) differs from one in [23] by a factor of [2"^^ — Ws4'^*i?i(l +a;s,0)]^. This 
extra factor agrees with the extra factor mentioned in the "iVote added^ of ref. [89] which 
also considered jet broadening in the context of SCET. Our result in eqn. (6.61) agrees with 
[89] up to the distinction between our oog given in eqn. (6.54) corresponding to factor rj in 
[89]. The distinction is conceptual and an important one when gauging scale dependence. We 
explain this below. 

In the formalism developed in [90], which was also used by Ref. [89], a single logarithm 
of appears in the combined result of the logarithm of IR sectors. The coefficient of 
the logarithm of is extracted and the logarithm is exponentiated when calculating the 
resummed cross-section. This coefficient is unique, however, the scale associated with in 
the logarithm is ambiguous in this formalism. Ref. [89] effectively make this choice same as 
the renormalization scale /i, hence ln(Q^//i'^) show up in rj, the factor corresponding to cOg- 
They could have made an alternate choice, for example ln(Q^/(Qi?r)^), since /i ~ QBt- The 
two choices differ in that they give significantly different estimate for /i- variation at NLL, 
one choice underestimates it while other overestimates. In our formalism such an ambiguity 
does not arise since u is unrelated to /x and can be varied independently to gauge the errors 
associated with the choice of this scale. We find that our combined /i and i^- variation errors are 
significantly smaller than the result of ref. [89]. This becomes important in phenomenological 
applications of the NLL distribution, for example fitting from the total jet broadening data. 

Ref. [89] state that our "analytic regulator" will not necessarily reproduce full QCD 
since we regulate in the effective theory. We have proven in this article that the regulator 
leads to a correct cancellation of ry divergences, preserves non-Abelian exponentiation, the 
soft-collinear gauge invariance and factorization in SCET. In SCET the analytic regulator 
leads to gauge dependent collinear and soft functions and breaks non-Abelian exponentiation, 
though it appears it can be modified to remedy these faults [69]. 

7 Application to Exclusive Processes and End Point Singularities 

We now discuss the application of RRG to the exclusive processes. There is a long standing 
problem in such processes, upon which, we hope to shed some light. In particular, many 
exclusive processes such as the pion form factor at large [91], the p — tt form factor [92], 
and the B — )• nlv form factor [93] are plagued by end-point divergences. In SCET these 
divergences arise in problems that fall with in the domain of SCETn. Thus to avoid this issue 
in a sensible fashion one is forced to remain in SCETj [94, 95] at the cost of loss of predictive 
power. If we could make sense of the end point singularities in SCETn, this would increase 
predictive power by allowing for one to write down the rates in terms of light-cone wave 
functions. Overcoming these end point singularities in a systematic fashion is thus highly 
desirable. 

^*Our regulator is not an analytic regulator by any definition that we are aware of, though it has the 
appearance of an analytic regulator at one loop for real radiation in the collinear sector only. 
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End point divergences arise when integrating over the momentum fractions in hght cone 
wave functions. Schematicahy these divergence ampHtudes are of the form 



where x is a momentum fraction (pn is a hght-cone wave function for hadron H, and C is 
a perturbative hard matching coefficient. It is often the case that C{x) is singular at the 
lower end-point. While we cannot calculate (pn from first principles, we can say something 
about its form in the limit where the hard scale is taken to infinity [96]. In particular when 
X approaches zero (pn vanishes linearly in x, whereas it is often the case that C{x) ~ l/x"^ 
leading to a divergence. Ref. [5] correctly account for this singularity as a rapidity divergence 
and performs a zero bin subtraction. However, in SCETn subtracting the double counting 
region can at best move the boundary between the soft and collinear sectors. The rapidity 
divergences and their accompanying logarithms still exist after the subtraction since they are 
associated with the boundary (see section (2)). As such, the discussion in [5] was incomplete 
in the sense that the method does not allow for a resummation. The crucial distinction, to 
be discussed below, between their regulator and the one employed here is that they imposed 
manifest boost invariance in each sector, whereas we explicitly break it with our regulator 
and only ask that the final answer be boost invariant. 

We follow [5, 21] in studying this issue within the context of a toy model of B decays. We 
will provide a proof in principle that exclusive i?-decays can be factorized in SCETn and all 
the logarithms can be resummed. We take the example of — )• iv'j where all fields are taken 
to be scalars. The physical case of fermions was considered in [21] but does not suffer from 
an endpoint singularity until the subleading order. The scalar case on the other hand does, 
and was originally thought to be non-factorizable due to these divergences [24] . Manohar and 
Stewart were able to regulate all the integrals and obtained the correct IR divergences after 
employing a zero-bin procedure. However, they were unable to resum the rapidity logarithms 
that appear in the ratio / fi~ , the scales associated with their rapidity regularization. This 
was due to a mismatch in hard logarithms and hard anomalous dimension. We will use the 
same set up as ref. [5] (see section VII A therein for details) but we employ the regulator 
introduced earlier in this paper. We will also include a region that was not discussed in the 
analysis of ref. [5] that we find essential for solving the problem in all frames of reference. 
We will ignore the wave function renormalization in our analysis as it is straightforward to 
include. 

At next-to leading order and in partonic approximation, the process B — )• ii^j is repre- 
sented by the Feynman diagram in fig. 11(a). In this process, observed in i?-meson rest frame, 
p is n-collinear with p~ tub and p2 = 0. While, £ is soft with ~ Aqcd and f = 0. After 
radiating the vector boson, a hard interaction takes place at vertex A and is only sensitive 
to the Lorentz invariant combination p^£^. The external on-shell particles (or particles with 
off-shellness at most Aq^j-,) are at vertices B and C. This full theory integral has IR diver- 
gences but is UV finite. To control IR divergences we choose to regulate the (scalar) quark 
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(a) B lu-f at LO 



(b) Sudakov form factor at LO 



Figure 11. (a) W^-boson and its decay is not shown. buW vertex is represented by •. u quark 
represented with the thick hne is the hard parton with ofF-shellness tosAqcd- Hard interaction takes 
place at vertex A and external on-shell partons are at vertex B and C. (b) Corresponding Sudakov 
form factor demonstrating hard interaction at vertex A. On-shell rt-coUinear parton is at vertex B and 
on-shell n-collinear parton is at vertex C. 

BC with a non-zero mass m, as was done it ref. [5]. We must take m ~ Aqcd to ensure proper 
scaling of the IR physics. This integral has a double logarithm, ln'^{p~£~^)/m?, that is large. 
A proper factorization should separate the hard scale p~i~^ from the non-perturbative scale 
m and provide a method for resummation. We will now show that this problem is identical 
to the Sudakov form factor at LO as far as the momentum flow is concerned 

We make a boost to a frame where £^ = p~ ~ \/fnB ^QCB- We will refer to this as 
the symmetric frame. In this frame, diagram in fig. 11(a) immediately maps to the familiar 
picture of the massive Sudakov form factor at one loop shown in fig. 11(b), where the hard 
interaction is at vertex A and (anti-)collinear on-shell quarks are at vertices B and C with 
p~ = = Q. The hard scale here is p^l'^ = and IR scale is set by the mass in of 
the vector boson exchanged between two quarks. The full theory integral is the same in the 
two situations when all particles are replaced by the scalars as in the toy example considered 
here. Previously we factorized the Sudakov form factor into (n)n-collinear and soft regions 
corresponding to the situation when the loop momentum k becomes (n)n-collinear or soft. 
The same factorization applies to -B-decays in the symmetric frame with the identification, 
~ ^^bAqcd, which then yields = m'^/Q^ ~ Aqcd/"i_b- 

Ref. [5], analyzed the problem in the S-meson rest frame (or lab frame), but did not 
include the region that corresponds to the soft region of the symmetric frame in their analysis. 
This region looks like a coUinear mode in the lab frame but with a lower rapidity compared 
to the coUinear modes that have the same scaling as the photon's momenta. In fig. 12, we 
show the three IR modes required for this problem in three very different frames of reference. 
To avoid confusion, we will refer to the modes as left, center and right modes corresponding 
to their location in the mode diagram. The three frames of reference considered are the 
symmetric frame, the lab frame in which i?-meson is at rest and the super-boosted frame 
where i3-meson itself is n-collinear. In the super-boosted frame both the initial parton and 

■^^A similar argument was used in [24], however these authors worked in SCETi power counting with off- 
shellness regulator which lead to the erroneous conclusion that the amplitude did not factor due to the existence 
of so-called messenger modes which are just the boosted ultra-soft modes. 
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Figure 12. Modes and Frames: relationship between symmetric, lab and super-boosted frame is 
shown along with the modes in each frame with their location on the hyperbola. 



the final photon are n-coUinear with diff'erent rapidity hierarchy. We will now demonstrate 
a factorization in the symmetric frame which essentially carries over to all other frames. 
The calculation of the operators changes between frames along with corresponding rapidity 
logarithms but the end result after the resummation stays the same. To get the same result 
in all frames it is crucial that our rapidity regulator be employed in a consistent fashion. 

For brevity, we will not present the details of the operators and factorization here and 
refer the reader to ref. [5]'^^. After suppressing an overall factor of ieg'^G/{p^ i^), we have 
the full theory integral for the diagram shown in fig. 11(a), 

^e^^^ I- ^ 1 
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where we see the large double logarithm that requires resummation. Note that there is no 
UV divergence in this full theory result, and that all the IR divergence are regulated by m?. 
In the symmetric frame, just like the Sudakov problem, the factorization is accounted for by 
the following integrals, each corresponding to a different region and operator, 
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(7.3) 



The diagrams corresponding to each of these integrals are shown in figs. 13(b), (c) and (d) 
respectively. They collectively provide renormalization to operator Oa shown in fig. 13(a). 
Here fr^^k) (shown in Table 1) corresponds to the additional rapidity regularization required 
to evaluate the integrals which arises from minimally regulating the corresponding operator 
by inserting an appropriate factor of fn therein. The presence of the regulator breaks the 



■^^The integral corresponding to the missing region and operator in this reference will be referred here as 
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(a) LO EFT operator 




4 



(c) k is soft 




— -<- 

(b) k is n-collinear 




) W\AA/ 



(d) k is n-collinear 



Figure 13. Diagrams in the effective theory contributing to B-decay in the symmetric frame. Each 
(g) corresponds to a different operator. 



boost-invariance in the otherwise invariant integrands. Therefore, in aU frames generated 
by boosts along n the same integrands arise but the regulator takes different forms in each 
frame. Our regulator is simply the correct limit of vP'v'^\'ik^\~'^ in each sector, where w is the 
bookkeeping parameter discussed earlier. In Table 1 we show the regulator in each frame for 
each mode. We will now discuss the calculation and renormalization in each frame. 

7.1 The symmetric frame 

Using the rapidity regulator as shown in Table 1 for the symmetric frame and set t(; = 1 we 
get, 
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symmetric frame 
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super-boosted frame 
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Table 1. /„ for each mode in each frame. 
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We immediately notice that the rapidity divergences cancel out in the sum of the three 
integrals, giving for the total bare effective theory contribution 
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The UV divergences in this result should be canceled by the counter term of the LO operator 
Oa- Therefore we obtain'^'' 
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for the counter term of Ja, the matching coefficient to Oa- The renormalized EFT contribution 
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Subtracting this from the full theory result of eqn. 
coefficient, 

eg'G 



IJL- 12 
(7.2), we get the one-loop matching 
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First thing we note is that the divergences in the counter term (7.6) are consistent with the 
logarithms in the matching coefficient J a, that is anomalous dimension obtained using eqn. 
(7.6) can be used to resum logarithms in Ja- Secondly, each operator corresponding to the 
left, right and center modes can be run independently in u to resum the rapidity logarithms 
in the infrared sector. Apart from the technicalities of operator mixing, the running strategy 
works the same as in the case of the Sudakov form factor discussed earlier. Therefore we 
have shown in principle that exclusive -B-decays can be factorized and resummed in SCETn, 
contrary to the previous claims [24]. 

7.2 The super-boosted frame 

Using the rapidity regulator for the super-boosted frame as shown in Table 1 we need to 
calculate /center and Jioft only, as /right is the same in all frames. 
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Note that in this toy example couplings e and g are quantities of mass dimension one. 
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where /center vanished because the last integral in was scaleless. Note that the left integral 
now has exactly the same ry-divergence as it was in the sum of /center and /left in the previous 
case, so are the structure of i/-logarithms. Therefore, the sum of three sectors, still 
yields the same result as in eqn. (7.5) of the symmetric frame. 



7.3 The lab frame 

Using the rapidity regulator as shown in Table 1 for the lab frame we note that /right is 
same as in eqn. (7.4) and /center = as was in the super-boosted frame. It only remains to 
calculate /left, whose exact evaluation is cumbersome. Therefore, for this frame we only give 
the ry-divergence structure of the left integral^^, 



/left 



|2fe 



{2ttY [A;2 - £+k- + iO] [A;2 - m2 + iO] [-p-k+ + iO] 



(7.10) 
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+ r/-finite 



Details of this calculation are shown in app. G. This is exactly the divergence required to 
cancel the rapidity divergence in /right- 



8 Conclusion 

In this paper we have presented a formalism which allows one to factorize and resum observ- 
ables which are sensitive to soft recoil. These observables fall within the confines of SCETn 
in which soft and coUinear modes have the same invariant mass scalings. It is because of 
this equality of scalings that one runs into rapidity divergences which force us to introduce 
a new regulator with an associated scale. We presented a proof that while individual sectors 
have rapidity divergences when one sums over sectors these divergences cancel as they must, 
since they are an artifact of factorization. The sectors contain soft and collinear function 
which are gauge invariant, and process independent. For transverse momentum distribution 
we are able to define a gauge invariant and universal transverse momentum dependent parton 
distribution function. Once the regulator is implemented one can sum the rapidity logarithms 
by use of the rapidity renormalization group, which corresponds to sliding the cut-off which 
separates collinear and soft modes on the mass shell hyperbola. 

We demonstrated our formalism by showing how one can sum the logarithms in the 
massive Sudakov form factor, as well as in the Higgs transverse momentum distribution and 
jet broadening. In the case of the Higgs distribution we give a generalized factorization 
theorem which goes beyond the classic CSS result in that it allows for jets in the central 
region. At leading order in the matching at the hard scale our results reduce to those of CSS. 
We also showed how our formalism can be used to renormalize exclusive processes with end 

■^^Naively, one may expect that this integral had an overlap with the central region. The overlap integral is 
obtained by taking k to be collinear (fc/x ~ Q(A^, 1, A)) in the integrand which exactly reproduces /center of this 
frame once the regulator is transformed as |2fc3|~'' — )■ Thus, there is no overlap between the regions. 



point singularities, which ahows one to complete the original calculations of Manohar and 
Stewart who showed how to sensibly handle the end point divergences using SCETn and the 
zero-bin subtraction method. 
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A Gauge Invariance and Rapidity Regulators 

In this section, we will prove the gauge invariance is not spoiled by regulating rapidity di- 
vergences concentrating on the generalized Wilson line regulator (rj). Generalizing the proof 
to the delta regulator follows in a simple fashion. We also show that regulating rapidity 
divergences with the rj or delta regulators is consistent with non-Abelian exponentiation. 

A.l Regularization at Higher Orders in the Sudakov Form Factor 

Given our operator definition, which we will see will have to be slightly amended when going 
to higher orders, it is not clear at all that gauge invariance is retained once the regulator is 
inserted. However, a simple argument shows that, at least within any set of covariant gauges, 
the regulator will preserve gauge invariance. Consider first the case of one gluon emission. 
In this case the gauge dependent piece of the propagator will generate an extra factor of 
n • k which then eliminates any rapidity divergence arising from that term. As such, for the 
integral involving the gauge piece we may set r/ to zero and gauge invariance follows. 

At higher orders we must modify the regulator in order to manifestly preserve both gauge 
invariance and eikonal exponentiation. We begin by considering the renormalization of the 
soft function which, as opposed to the jet function, is a pure Wilson line. Let us recall some 
basic facts about Wilson lines and their renormalization: 1) The anomalous dimensions of the 
cusped Wilson line is at most linear in logarithms [37, 38]. 2) The result exponentiates at the 
level of the integrands [34-36], with each color Casimir appearing only once. At each order 
there are a set of graphs which are two eikonal line irreducible (2EPI), i.e. they can not be 
disconnected by cutting two eikonal lines. These graphs generate a Casimir which does not 
appear in any lower order graphs. The sum of these graphs is called a CWEB. Thus in the 
exponent only CWEBs appear. These two facts imply that the sum of integrands which form 
a CWEB have no higher order power beyond 1 /e^ for UV or 1 /rj for rapidity. Poles of order 
l/?7^ would violate condition 1) since all integrals have UV divergent transverse momentum 
integrals that are regulated by dim. reg. thus leading to anomalous dimensions which are 
not monomials in logarithms. 
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Figure 14. 



To prove gauge invariance we should choose our regulator so that we may sum the di- 
agrams which are part of a CWEB. Doing so allows us to take advantage of the lack of 
sub-divergences in this subset of diagrams. Thus instead of regulating the individual mo- 
menta, we will regulate the group momenta n- kg. Where kg is the total momentum flowing 
into a CWEB on the eikonal line . So we will re- write our regulated Wilson line as 



Wn= ^ exp 



perms 



n ■ Ar. 



Sn= exp 



gw \ n ■Vg\ ^ 
~n- V 

gw \2Vsg\-^/\ 



(A.l) 



(A.2) 



perms 

Note that non 2EPI diagrams will thus be regulated differently then 2EPI diagrams. Consider 
the 2EPI diagrams which contribute to the O(a^) CWEB shown in Fig. 14. The sum of 
the integrands contains no sub-divergences. By regulating only the group momentum the 
integrands may still be added. So the sum of the integrands have the general form 



ICWEB^ j d%d'^kD{k',klk-kg)-^ 



\k^g\ N{n ■ k,n ■ k,n ■ kg,n ■ kg) 

{h ■ kg + ie){n • kg — ie) D{k'^, k"^, k ■ kg,n ■ k,n ■ k,n ■ kg, n • kg) 

(A.3) 

Here we have explicitly pulled out the first eikonal propagator. Note that in all contributions 
to the CWEB the final gluon attaching to the eikonal line will carry the same momentum as 
the eikonal line itself. Therefore, if we now consider a general covariant gauge with polariza- 
tion sum 

Y.',e.-9,u + {l-0'^- (A.4) 

The gauge dependent piece will necessarily cancel the final eikonal propagator and since the 
CWEB has at most an order I/77 rapidity divergence the gauge dependent piece will always 
be finite. Given that the gauge dependent pieces are finite, we may set r] to zero in those 
contributions, leading to a gauge invariant result. 

Now let us consider the set of diagrams which are not 2EPI such as the one shown in 
Fig. 15. The order CWEB involves two additional sub-diagrams corresponding to the 
vacuum polarization and the Y-graph as in Fig. 14 as well. Eikonal exponentiation implies 
that the sum of this diagrams factorize into a product. That is, this contribution comes 
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Figure 15. An example of a diagram composed of two CWEBS. The group momentum for tlie 
CWEBs is labelled kgo- 

from interference terms in the expansion of the exponent and will arise with a color factor 
CpCA- The rapidity regulator preserves this property by regulating the group momentum of 
the individual CWEB independently. That is the integrand for this diagrams will be of the 
form 



By defining the regulator in this way we preserve the eikonal identities utilized to show that 
the sum of graphs yields the product of a one and two loop integral. 

The renormalization and gauge invariant nature of the jet functions follows by a similar 
line of reasoning discussed above. However, there is a crucial difference between the jet and 
soft functions. Since the jet only involves one Wilson line the notion of a 2EPI diagram is 
no longer applicable. However, as will be discussed in the next section, it is simple to define 
a CWEB for the jet function when noting that every CWEB diagram in the soft sector has 
a parent diagram in the full theory which has a collinear limit. There is also the matter of 
pure self energy diagrams. Such diagrams on the non- Wilson line part of the jet function, 
obviously have no rapidity divergence since there are no eikonal lines involved, i.e. they are 
identical to pure QCD. All of the eikonal self-energy graphs vanish after including the proper 
soft-bin subtraction, since the full theory self energy graphs are reproduced by the purely 
collinear diagrams. This can be seen from explicit calculation as well. 

A. 2 Regularization of Generalized Soft and Collinear Functions 

For generalized observables the soft functions are more complex then the one which arises in 
the Sudakov form factor. Typically we are interested in soft matrix elements which arise from 
the amplitude squared and where we measure some aspect of the soft radiation that can be 
written as 



where p stands for a set of momenta components which scale as some non-zero power of A. 
The soft function will necessarily be accompanied by at least one collinear function which can 
be written as 




(A.5) 



S{p)^{0\SnSlS{p-V)SlSn\0) 
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(A.7) 
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The operator O here is either a quark or a gluon field and Wn here is a cohinear Wilson line 
in either adjoint or fundamental representation depending upon the case. The delta function 
measures the kinematic quantity of interest and ip is either a vacuum ( "Jet Functions" ) or a 
hadronic state ("Beam functions", "TMDPDF's" or generalizations thereof). 

All of these functions will contain rapidity divergences in SCETn upon factorization. It 
is only the combination that must be (rapidity) finite. The finiteness of the total result is 
predicated on the fact that the diagrams are regulated in a consistent fashion. In particular, 
as touched upon in the previous section, the regulated EFT diagrams must come as pieces 
of the asymptotic expansion of the "regulated" full theory diagram. We use quotes here to 
remind the reader that the full theory has no rapidity divergences. It should be clear that 
there are many ways in which to insert our rapidity regulator into the full theory diagram. 
However, regulating the full theory diagram does not guarantee that the EFT diagrams will 
be regulated. One can imagine that the full theory, upon expansion, has rapidity divergences 
in sub-diagrams and that a particular choice of regulator could in general lead to unregulated 
EFT diagrams. Furthermore, we should choose a method such that our generalized function 
are universal. That is, they should be process independent, though they will always be scheme 
dependent just as with any parton distribution function. Finally, we need to ensure that the 
regulator preserves eikonal exponentiation. These criteria are not logically independent so it 
should not surprise the reader that its relatively simple to ensure that all are satisfied. In 
particular once we ensure that we preserve exponentiation, the other criteria are automatically 
satisfied. 

To discern a proper prescription we begin with the soft function where the notion of a 
CWEB is clear. As in the case of the Sudakov form factor we regulate the total momentum 
emitted in a CWEB. We will work in the Feynman gauge and the proof of gauge invariance 
(covariant gauges) follows by an the identical argument given for the Sudakov form factor. 
Note that for the soft function there is no need to ever calculate any diagram which is 
not a CWEB. At a fixed order a non-CWEB diagram will not contribute to the anomalous 
dimensions. Furthermore, for the purposes of matching, all non-CWEB diagrams can be 
determined by expanding the exponent, i.e. non CWEB diagrams are simply products of 
CWEBS (recall each CWEB is affihated with a color Casimir). Summing all the CWEB 
integrands ensures that the diagrams are marginally divergent. Thus, we regulate the soft 
function by inserting a factor of | k-^g \ ~^ where k^g is the total momentum flowing from the n 
to n side. The marginal nature of the divergence assures us that this choice of regularization 
is sufficient. 

We may associate a collinear CWEB with each soft CWEB. This is true despite the fact 
that there is only one eikonal line in the jet function. While there is no soft contribution 
to the collinear function, after soft-bin subtraction, we may consider the soft limit of the 
diagram for the express purpose of determining whether a diagram is part of a CWEB. This 
is exactly what happens in the soft limit of the parton distribution function [97]. Thus given 
our choice of regulator for the soft function we regulate the collinear function by inserting 
a factor of | kg± where kg is the total momentum fiowing off the Wilson line. Note 




Figure 16. The correspondence between the fuU theory diagram and the diagrams in the effective 
theory. The first column on the right is the soft function contribution while the next two columns 
correspond to the collinear functions in the n and n directions respectively. Not shown are mirror 
diagrams in which n n. The blob corresponds to possible hard pieces or soft non-hadronic (e.g. the 
Higgs) final states. The darkened lines are eikonalized. The incoming lines in the collinear functions 
are either quarks or gluons. 

that the rapidity divergences in the collinear function will exponentiate (as required by the 
cancellation of rapidity divergences), but the entire function, including UV divergences, does 
not. 

Let us illustrate how this works at two loops. On the left hand side of Fig. 16 we have 
a full theory diagram. We will ignore hard gluons since they are handled trivially. The full 
theory diagrams can be expanded around three effective theory diagrams (plus their mirrors). 
The possibilities correspond to two softs, one soft one collinear and two collinears (in each 
direction). The resulting factorized EFT diagrams are shown on the right hand side of the 
figure. The purely soft diagram corresponds to a contribution to the 0{a'^) CWEB. Note that 
at two loops the purely soft (or collinear) diagrams will not be marginally divergent, only the 
sum of CWEB diagrams will be. The sum of the diagrams has no rapidity divergence. Note 
that the factorization which occurs in the middle diagram in general will only occur once all 
diagrams have been summed over. In Fig. 17 we have a case where the collinear function is 
not purely eikonal. i.e. the gluon running straight across the cut is not part of any CWEB 
and multiplies the exponentiated CWEB. 

In the end the regularization is straight forward. In the soft function insert a factor of 
I 2A;3 I"'' where ^3 is the total momentum in the spatially longitudinal direction flowing from 
n to n. In the collinear functions (whether they are beam or jet) insert a factor of | k± |~'', 
where k is the total momentum flow out of the Wilson line. This defines process independent 
collinear and soft functions, including TMDPDF's. Once any such function is calculated it 
can be used in any process as long as all the contributing function are calculated in the same 
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Figure 17. The gluon running straight across the cut is not part of any CWEB and does not 
exponentiate. 

scheme. Note that not all schemes will share this property of universality. For instance, 
suppose we chose to regulate the momentum flowing across the cut, in which case purely 
virtual diagrams would be set to zero. This scheme is attractive for its relative simplicity of 
definition. However, while this would give a sensible beam function or any such generalized 
PDF, it would not regulate any collinear function which corresponds to a vacuum matrix 
element since the total momentum across the cut will sum to a constant. 

Finally we must consider the case of self energy diagrams. As we argued at the end of 
the section on the Sudakov form factor, the self energies must be purely collinear since these 
graphs are identical to the full theory graphs. All eikonal self energies must therefore vanish 
after soft bin subtraction. 

A. 3 The Use of Other Regulators 

As was briefly mentioned previously, it is certainly possible to use other regulators in this 
formalism. We emphasize that the rapidity renormalization and resummation itself is reg- 
ulator independent formalism, just as the renormalization and resummation of traditional 
ultra-violet divergences are independent of dimensional regularization, but are often most 
conveniently performed in dimensional regularization. To be able to renormalize and resum 
rapidity divergences is also in the same sense a regulator independent procedure, though to 
implement in any given instance some regulator must be used. Given that, there are certain 
properties a regulator should fulfill: 

• Gauge Invariance 

• Preserve Non-Abelian Exponentiation 

• Have a universal definition for generalized soft and jet functions 

• Clearly delineates sectors 

As an example of a regulator that can be engineered to satisfy all four conditions, we 
consider the J-regulator of [32]. At first, the J-regulator does not manifestly preserve eikonal 
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exponentiation, as can be seen even at the Abelian level where the famous eikonal identity 
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a=l 
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is no longer obeyed due to the 0{6) terms. Nonetheless by adding all the elements of a 
CWEB only the 5 divergences of a single log will appear, and these 0{6) terms can be safely 
set to zero. That is to say, only the finite and log divergent pieces survive the 5 — )• limit: 
in a CWEB, higher order terms are never needed. Further, the arguments outlined above 
illustrating the rapidity integration finiteness of gauge dependent terms work as well with the 
(5-regulator: the regulator can safely be set to zero in such terms. 

The subtlety of the 6 regulator is in the jet sectors. Before one can renormalize the 6 
divergence in the jet sector, one must perform a "soft"-bin subtraction. This subtraction 
does not remove the rapidity divergence, but rather forces that sector to have the correct 
rapidity divergence [32]. This shifting of rapidity divergences in SCETn is analogous to zero- 
bin subtractions in SCETj. The zero-bin in SCETj enforces each sector to have the correct 
ultra-violet divergences, by removing the ultrasoft divergences which should not be attributed 
to that sector. Any rapidity divergence completely disappears (sector by sector) in SCETj 
(with the inclusion of the zero-bin) but ultra-violet divergences remain. This accentuates the 
fundamental distinction between the two theories. Whereas SCETi has mode factorization 
in the invariant mass, SCETn has mode factorization in rapidity. The evidence of such 
factorization in a perturbative calculation is the divergences arising in the integrals of each 
sector, whose renormalization allows for the resummation of large logarithms. 

After the soft-bin subtraction of the jet sector has been performed, one can renormalize 
both the rapidity and ultra-violet divergences. In the case of the 5-regulator, this will trade 
the regulator 5 for an arbitrary parameter analogous to our scale i'. The regulator itself 
should be formally removed, and and arbitrary scale corresponding to the renormalization 
point introduced. At this point one would see that the inclusion of the terms with a positive 
power of 5 in (A.8) was unnecessary. 

B Soft-Bin Subtractions 

In SCET the zero bin subtraction plays an important role even when it formally vanishes. 
In SCETi collinear modes have invariant masses which are parametrically larger then the 
ultra-soft modes. Thus when we calculate a collinear loop we should not expect an ultra-soft 
divergence, only a collinear one. Invariably when calculating loops of the collinear modes one 
does generate Ultra-soft (non-collinear) divergences which should be attributed to the ultra- 
soft sector. By doing zero-bin subtraction (by zero-bin subtraction, we mean the subtraction 
of the ultra-soft contribution to an integral, where the loop momenta is taken to scale as 
Q(A^, A^, A^) in light-cone coordinates) this IR divergence is morphed into a UV divergences. 
This is sometimes called the "pull-up" mechanism [5, 98]. 
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In SCETn since the soft and collinear modes are on the same mass sheh hyperbola we do 
a soft-bin subtraction. In doing a "soft-bin" subtraction we subtract the soft region, where 
the loop momenta is taken to scale as Q{X, X, X). The soft-bin of the collinear mode serves 
the purpose of cutting off the integral at the proper spot on the rapidity hyperbola. This is 
most clearly seen when one regulates with a d regulator as was shown in [32] . In this work the 
authors introduce a different regulator for each mode. The particle i propagator get replaced 
with 



1 1 



(B.l) 



{pi + ky - mf {pi + ky - mf + Aj ' 

For the two collinear modes, we will call the momenta p±. A is designed to distinguish (cut- 
off) the + momenta from the soft region. However, the emission of a kj^ gluon of the p- lines 
leads to an eikonal propagator P_ of the form 



1 

k+ + 



V- 



These emission build up a Wilson line in the + sector, with an apparent violation of factor- 
ization, since the collinear + sector now depends upon the wrong cut-off (i.e. the it sector is 
sensitive to A/p_). This cut-off must be pulled- up to the correct sector by the soft-bins and 
this is exactly what the soft-bin accomplishes [32] ^'^ . 

When using an r/ regulator the soft-bin is scaleless and is thus vanishing. Its formal 
subtraction from the collinear integrals has the same physical effect as in the case of the delta 
regulator. To see this more clearly we may work with both the r] and 6 regulators where the 
limit (5 — )• is take first. In this case the soft-bin subtraction acts to eliminate the "wrong" 
delta cut-off and replace it with the RRG scale v. 

Finally, it is important to emphasize that the collinear matrix elements and Lagrangian 
are well-defined only with the implicit soft or zero-bin subtractions [5, 32]. In particular, the 
subtractions are necessary to make the matrix element gauge independent [32] in each sector. 
This role of the soft-bin is unrelated to the issue of rapidity divergences, as we have proven 
the rapidity divergences to be gauge independent. These soft-bin subtractions can often be 
related to a matrix element of soft Wilson lines, and often these matrix elements of Wilson 
lines are the inverse of the soft-function found in factorization. In this guise they appear in the 
modern CSS formalism [70], where the inverse soft factors play a similar role to ensure gauge 
independence. However, as is the case in section 6, the inverse of the soft-function cannot 
be always identified with the soft-bin subtraction (even if the subtraction is representable as 
some matrix element of Wilson lines) , due to the differing kinematics and phase spaces of the 
soft function and jet functions. 



^^Note that there is no coUinear subtraction from the soft because the soft diagrams are insensitive to the 
large scales and there can be no overlap. This is not true if there are external particles with soft momenta as 
in the B meson decays discussed in [5] and in section (7). 
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C Higgs b calculation 



C.0.1 6-space TMDPDF 

For completeness we note the impact parameter space results, where the transform is strictly 
two-dimensional. At tree level, we have: 
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(C.l) 



At one loop, we get: 
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Ignoring the traceless term and expanding in r] and e give: 
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We can find the appropriate renormalization factor through the matching procedure to the 
PDF (see (5.49)): 
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The evolution equations satisfied are 
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CO. 2 6-space Soft Function 

In impact parameter space we have for the tree-level result and one-loop correction: 

S^^\b) = {2Trf (C.IO) 
= i2nf '~"~J^'' e-^^.V^-fe^-H-. ^(-^ -|) ^^[^^^"^^ (C.ll) 

Separating out the divergent and finite parts gives: 

(C.13) 



4^(^1 = 1- — 
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The evolution equations satisfied are 
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At one loop we have 
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iy-^S{b,ii,iy) =^srap{bij)S{b,ii,v) (C.16) 
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D Jet Broadening Resummation in Fourier-Laplace Space 

For completeness, we will also perform RRG in Fourier-Laplace space, as in comparison to 
some older literatures. The factorization theorem in b-r space can be written as 
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The soft function in 6 — r space is given by 
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The renormahzed soft function to NLO can be obtained by taking the finite part of Eq.D.3 
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D.l Renormalization Group Equations 

From the divergent part of theEq.D.3, we can get the renormahzation factor 
Zs(rR,TL,bi,b2) = 1 + „ n n H — — In 
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The /U-anomalous dimensions follow in standard fashion, and the i^-anomalous dimensions can 
be obtained by using Eqs.4.19 and 4.26 
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To eliminate the large logarithms we may run in both // and v to some fixed scale. That is 
we may write 



S{n,v) = Usifj-,fj-s;j^s)Vsi'^,i's;tJ')Sifj-s,'^s) 



(D.8) 
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where 

Vs{i^,i^s; ^J') =Exp 
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The hard function, with no running in nu, has the standard evolution, 

H{Q,^l) = u{^^,^XH)H{Q,^XH): 



UHifJ',fJ'H) = Exp < 2 
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D.2 Resummation at NLL 



In 



2/3o Vas{^J^H) 



Resummed broadening distribution to NLL in the Laplace space including both tradition RG 
and RRG depends on tree- level hard, jet, and soft function, plus two-loop cusp and one- loop 
non-cusp hard anomalous dimension, together with one-loop soft z^-anomalous dimension. 
The hard function, with no rapidity divergence and therefore no u running, are as standard 
hard function. 

The resummed left and right broadening spectrum to NLL can be written accordingly as 

1 dV^LL 
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in which the jet function in 6 — r space as defined in Eq.D.2 calculated at tree level as 
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We can then write Eq.D.15 as 
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The left and right broadening distribution can then be calculated by performing inverse 
Laplace transformation 
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For total jet broadening, we have 
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To compare with standard the Total Jet Broadening {Bt) definition[3], where Bt 
write the broadening distribution as 
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aodBr " T{2ujs) Bt^BtQJ 2-^'^- \ l + Us 
equivalent to that is derived in the physical px-e space using the relationship between Hy- 
pergeometric and incomplete Beta function. 

E Transforms 

We collect here some transforms needed for resummations or deriving plus distribution prop- 
erties for both jet broadening and higgs. 



E.0.1 Fourier Transforms For Higgs spectrum 
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Expanding both sides in a and identifying powers of a gives the transforms of the plus 
distributions L„(/Lt,px). 

E.0.2 Transforms For Jet Broadening 

Transforms: 

dxx"J^(x6)e--- = 2-^b^r-"-^-^ ^^^ + ^ + ^^ 2i^i Q(l + a + ^(2 + a + ^); 1 + /?; 

(E.3) 

Inverse Laplace Transforms: 



Inverse Bessel Transform: 



JO 



r(w) (e2-p2)] 



F Plus-Distributions over Vector Domains 

Both the jet broadening and Higgs transverse momentum cross-sections involve convolutions 
in vector quantities. Expressing the renormalized functions directly in momentum space 
requires plus distributions that automatically perform the subtractions necessary to render 
convolutions finite. We give definitions for these plus distributions, to make transparent the 
renormalization scale dependence of the functions, and give useful identities for their manip- 
ulation. For the purposes of this paper, we consider mostly one class of these distributions: 



(F.l) 



Intuitively, these are distributions that render integrals over a vector domain convergent when 
weighted with a well-behaved function. This is accomplished by subtracting away an integral 
over a disc of radius X/i (formally, P^/i = {p '■ \p\ < X/j.}) about the origin in the p-space. 
Thus: 



(27r)2' 



C^U,q;\)=0. (F.2) 



From this class, there are important limits that appear frequently, and are related by deriva- 
tives. First we define the notation 

£-(m,p;A) =£?(^,p;A), (F.3) 

Cn(^IJ,,p;\^ = Cn(^iJ-,p;\^ . 
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Then we have the relation 



hm — — - 



(F.4) 



Henceforth, we wih also assume the notation that when the argument A is absent it implies 
it is set equal to 1. 

In what follows, we will give a definition based on dimensional regularization since it is 
extremely useful to perform multidimensional integrals. Any other intermediate regularization 
is also fine, since all divergent behavior cancels in integrals with well-behaved test functions, 
and so any such regulator will cancel after performing the integral. We have checked, for 
instance, that the use of a limit based definition with a mass regulator leads to the same 
expressions. 



F.l Definition in Dimensional Regularization 

One simple way to define a plus-distribution with vector arguments is dimensional regular- 
ization. This is especially appealing given the well-developed nature of the technology and 
its familiarity. Assume g{k) is a reasonable test function, that is, has no singularities at 

the origin, and falls off at infinity fast enough. For physical applications, these criteria are 

I A 



easily met. The plus-distribution, 



for a function / that has at most a power-like 



singularity at the origin and the boundary condition on the disc V\^^ is defined as: 
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We use the notation: 



A] = /Li 



-2e 
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^2+2. 
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B,[f;f,]=B,[f;f,,l]. 



(F.6) 
(F.7) 



Formally, all integrations are in two dimensions as long as the [•]+ symbols are used, though 
implicitly the integrations are in 2 + 2e dimensions where analytical continuation on e is 
assumed as always. /^"^"^ that appears on right hand side of eqn. (F.5) is present only 
to formally control the dimensions'**'. The use of makes the transition between these 
measures less cumbersome. Put simply, = (Itt)"^ 6^'^\p) or {2tt)'^^'^''6^'^^'^'^\p) depending 
on the context. Further, restriction to discs about the origin is not necessary, but merely 
convenient. Any simply connected region containing the origin can serve to define the plus- 
distribution. Finally, the boundary term, B^, by definition gives the following identity: 



Xjj, 



(2vr)2 



0. 



{¥1 



Any other momentum scale will equivalently serve the purpose. 



- 74 - 



when we replace the test function g{q) = Vx^xiq)- Two different boundary conditions for the 
plus-distributions are related by: 



f{fJ',P 



lJB,[f;fi,X,]-B,[f;i^,X2]) 



(F.9) 

For the distributions in (F.3), we can give explicit expressions for the boundary terms in 
dimensional regularization^^ : 
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where on the right-hand side, e appears explicitly and p is understood to be a vector in 
2 -|- 2e dimensions. It is important to emphasize that the limit e — )• can only be taken in 
expressions where the limit is manifestly finite. In most applications this occurs only after 
evaluating integrals like (F.5), or for instance in the difference between two different boundary 
terms as in (F.9). 

The power-law distribution A) becomes especially simple with the A = oo bound- 

ary condition. In this case, the boundary term vanishes, and we have the identification 
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2n l+c 



1 //i2xl+ 



/i2 \p2 



It is useful to give the expansion of this distribution as a power series in a 
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27r/Li2 
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(F.15) 



F.2 Rescaling 

We exhibit the scaling identities analogous to those defined for scalar domain distributions, 
for example see [99]. For the distributions and Cn{jJi-,^, these are essentially the 

same as in the scalar case and take the following form in our notation: 



P' 



2a 



1 



27r 47ra 



n 

where nCm = r(l+m,jr(l+n-m) • 



m=0 



% ln"+^(p^) 
27r 47r(n + 1) 



(F.16) 
(F.17) 



'^Here we have pulled an explicit 1/2-k compared to definition (F.6) for convenience. 
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F.3 Derivatives and Integrals 

A few important derivative identities are: 



^£„(/i,p) = n£„,_i(^,p) . (F.19) 



2 d 

We can easily integrate the distribution C°^[^^p) over a disc of radius A/i by using the 
relations (F.8) and (F.9): 



^ A-2"-l) (F.21) 
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Similarly: 



F.4 Convolutions 

We have for the convolution of £"( and \fi,p 



"^'P Cjf,,p) = , 1/^" , ln"+^(A^). (F.22) 



/ 



+ {c/(a,/3)i?[£"+'^;/i] - B[£''; ^,]B[£^■, ^,]] (F.23) 

m r(l + « + /3) r(-«)r(-/3) 

Using the fact that logarithmic distributions are related to power-law ones via derivatives, one 
can use this convolution formula (F.23) to derive identities for the convolution of logarithmic 
distributions. However, care must be taken in the limit of vanishing power parameters. 

G Structure of divergence for an integral in sec. 7.3 

We wish to find the divergence of the following integral from sec. 7.3, 

^eJ2^ |2fc3| 

{AttY J {27ry [k^ - e+k- + iO] [A;2 - m2 + iO] [-p-k+ + iO] ' ^ ' 

Consider the auxiliary integral 

^(aux) ^ a fJ- n i " ^ l^'^sl ,p 2^ 

^^^^ (47r)^ y (2^)'^[-£+A;- + iO][fc2-m2 + iO][-p-A;+ + iO] • ^ ' 
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For the difference, /left — I\gf^^\ with some effort it can be shown that 



^ _ .(aux) _ -im^ e^'^g/i^V^ 



dks dkf (G.3) 



[k1 + m2) + A;,^ + m? ^kj + k'f + - k3 - ^ - iO 



+ r/-finite 



where the integral shown is the only integral that has a rapidity divergence. The divergence 
appears when ^3—7-00. So, as far as the rapidity divergences are concerned, it has the same 
asymptotic properties as the integral 



167r 2^? 7o + m?) y^kj + k'^ + m? 



(G.4) 

In other words, /left — Iieft^^ ~ -^asymp does not have any rapidity divergence. Since, ^ig^^^^ is 
the same as /center of eqn. (7.4) we conclude 
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